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Abstract 

In the large N c limit cold dense nuclear matter must be in a lattice phase. This applies also to 
holographic models of hadron physics. In a class of such models, like the generalized Sakai-Sugimoto 
model, baryons take the form of instantons of the effective flavor gauge theory that resides on probe flavor 
branes. In this paper we study the phase structure of baryonic crystals by analyzing discrete periodic 
configurations of such instantons. We find that instanton configurations exhibit a series of "popcorn" 
transitions upon increasing the density. Through these transitions normal (3D) lattices expand into 
the transverse dimension, eventually becoming a higher dimensional (4D) multi-layer lattice at large 
densities. 

We consider 3D lattices of zero size instantons as well as ID periodic chains of finite size instantons, 
which serve as toy models of the full holographic systems. In particular, for the finite-size case we 
determine solutions of the corresponding ADHM equations for both a straight chain and for a 2D zigzag 
configuration where instantons pop up into the holographic dimension. At low density the system takes 
the form of an "abelian anti-ferromagnetic" straight periodic chain. Above a critical density there is a 
second order phase transition into a zigzag structure. An even higher density yields a rich phase space 
characterized by the formation of multi-layer zigzag structures. The finite size of the lattices in the 
transverse dimension is a signal of an emerging Fermi sea of quarks. We thus propose that the popcorn 
transitions indicate the onset of the "quarkyonic" phase of the cold dense nuclear matter. 
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1 Introduction 



The phase diagram of QCD at finite temperature and chemical potential has a very rich and non-trivial 
structure. It is believed to possess a number of interesting phases and transitions between them (for reviews 
see [TJ [2l [3j [4j [5] . ) At high temperature QCD must exist in a deconfined phase with restored chiral symmetry, 
while at small temperature and high density it is believed to be in a color superconducting phase. These 
two regimes attract a great deal of attention because they are "accessible" experimentally in the heavy ion 
collisions and interior of compact stars0 Nevertheless the most challenging and perhaps interesting and rich 
physics resides outside these limits. 

Despite decades of study of QCD and nuclear physics, the current theoretical understanding of the QCD 
phase diagram is still very limited. Due to strong coupling the direct approach to QCD is only possible 
via lattice methods. Strictly speaking the latter are only applicable in the zero density case, since finite 
density triggers the notorious sign problem. Fortunately some technologies were developed recently in the 
lattice community, which allowed to extrapolate the zero density results somewhat off the zero density 
regime, though not too far. It turns out, however, that there is another loose end one can follow to untie 
the tangle: the large density case QCD is expected to become perturbative due to the presence of a large 
scale - chemical potential. The analysis of this phase can then follow the one of weakly coupled BCS 
superconductors. Accidentally this two theoretical handles we have to study the phase space are applicable 
precisely in the accessible experimental regimes. Such a situation looks in fact unfortunate because it leaves 
us no efficient tools to understand the intermediate regime, for which we have to resort to all kinds of 
bottom-up or effective approaches. 

In this work we will be interested in the properties of nuclear matter, which lies precisely in the inter- 
mediate part of the phase diagram. Therefore we have to choose an effective tool to study it. One of the 
relatively new effective and efficient tools was brought up by the ideas of holographic correspondence, or sim- 
ply holography, which received a precise formulation with the advent of AdS/CFT [BJ. Holography provided 
a new perspective for studying and understanding gauge theories in the regime of strong coupling. Among 
its many interesting developments it contains encapsulated toolkits to study theories which are cousins of 
QCD and hadron theories. 

In this paper we have chosen to work with a generalization of Sakai-Sugimoto model [7], which is so 
far the best-known holographic model of hadron physics. It is engineered from the following blocks. To 
construct the pure gauge sector one uses Witten's model of confining non-supersymmetric gauge theories [8] . 
Namely one works with a near-horizon limit of the theory on the D4-branes compactified on the circle 
in type IIA theory. The flavor sector is added following the ideas of Karch and Katz [9], by introducing 
additional D8-branes in the geometry of the D4. Notice that flavor branes are typically added assuming 
quenched approximation Nf <C N c , to neglect the backreaction on the background geometry. Sakai and 
Sugimoto were first to notice that D8 and anti-D8-branes smoothly connected in a U shape manner, model 
the phenomenon of chiral symmetry breaking in holography. Following a pioneering work by Kruczenski 
et al. [TO] , mesons can be analyzed as fluctuations of the fields in the effective low-energy theory of the 
flavor branes. Various achievements of holography applied to meson physics are reviewed in [11] , In turn 
baryons are introduced as other D4 branes wrapped on compact cycles according to the recipe suggested 
by Witten |12j . Since the work of Hata et al. |13j baryons in Sakai-Sugimoto model were studied in many 

1 Here the quotes stress our poor control and understanding of the experimental data. 
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works. See [TTJ and a more recent review []3] for a (non-exhaustive) list of references. 

The study of QCD phase diagram has recently drawn a revived attention mainly due to experiments on 
heavy-ion collisions and progress in theoretical understanding of QCD phases. A number of papers have 
appeared that studied the phase diagram through holography including [15] [16j EH OS US] ■ In holography 
though for the sake of feasibility one is bound to study a specific limit of gauge theories, namely large number 
of colors N c and large 't Hooft coupling A. In this limit the phase diagram gets important modifications. 
For the nuclear matter phase, which is of main interest for us, this means two major changes. First, in 
the large N c limit the phase of nuclear matter is a crystal, as opposed to ordinary nuclear matter. This 
happens because the potential energy of baryon interaction scales as N c , while the kinetic energy is only 
N^ 1 . Second, in the large A limit baryon mass tub scales as A, but baryon interactions only as A . As a 
result baryon density grows very fast as a function of the chemical potential above the critical value /i c ~ ms- 
At densities such that baryons start to overlap, quarks no longer know which baryon they belong to and 
the matter essentially becomes a quark matter. Therefore crystalline nuclear matter occupies only a narrow 
region of the phase diagram plotted in units of /x c and it disappears in the precise A — > oo limit. In the case 
of Sakai-Sugimoto model one can find a A — > oo diagram in the work of Bergman, Lifschytz and Lippert [17] . 

We are going to zoom on the phase diagram of Bergman et al. around /i = /i c and T = in order to 
study crystals of holographic baryons. A general idea of construction of baryonic lattices in Sakai-Sugimoto 
model was already outlined in the work of Rho, Sin and Zahed |20j in relation to skyrmion lattice, which we 
discuss below. We are going to put the ideas of Rho et al. on a solid ground by working out several examples 
of lattices explicitly. In Sakai-Sugimoto model baryons can be identified with instantons of SU(Nf) gauge 
fields that live on the world volume of flavor brancs. Therefore a description of a baryon lattice requires 
a construction of the corresponding instanton solution of the gauge field equations of motion. Notice that 
these are curved space instantons, so the full solution for them is not known. Instead an approximation is 
derived employing first few orders in 1/A expansion, following 13 . 

In the first order in 1/A expansion the curvature can be neglected. The baryon is a flat space-time 
BPST instanton. Such a solution is characterized by moduli, such as the instanton positions, radii and 
orientation. To compute interaction energy of instantons one has to go to the next order in perturbation 
theory. When interactions are turned on the moduli are partially stabilized. In particular the instanton 
acquires an equilibrium size. One can then quantize the resulting Hamiltonian to derive the baryon spectrum. 
Thus to find the ground state of a baryon system one has to find the relevant instanton solution, compute 
corrections to the BPS energy due to curvature and self-interactions on the BPS solution itself and minimize 
the net energy with respect to the moduli. 

We will follow this general strategy to address the following question. The energy of a single baryon is 
minimized, when the baryon is sitting in the lowest point of the gravity potential (bottom of the U shape 
profile) . When a finite density baryon configuration is considered the repulsive interaction between baryons 
should in principle be able to push them out from the bottom of the gravity well. In other words, suppose 
we consider a 3D lattice of instantons and squeeze it to an arbitrarily large density. Will the instantons stay 
bound to 3D or will they pop out from their 3D alignment? 

There are several pieces of motivation for posing this question. First, this question appears to be natural 
from the point of view of holography. If the above expectations are confirmed then the next question is of 
course the interpretation of this result. On one hand the answer to these questions is partially known thanks 
to the work of Rozali, Shieh, Van Raamsdonk and Wu |18j . They have shown that a uniform distribution 
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of instantons in 3D will necessarily incur a mesoscopic instanton density in holographic dimension as well. 
Finite instanton distribution must then be interpreted as a Fermi sea (of quarks, since the baryons in this 
picture are classical, neither fermions, nor bosons) and its edge as a Fermi surface. On the other hand the case 
of a uniform density is slightly different. A uniform density implicitly assumes the situation of overlapping 
baryons, when the quarks can no longer distinguish the baryons and start behaving in a qualitatively new 
way. Hence it will be interesting to understand physics in the non-uniform case as well. It turns out that 
the latter case has a rich structure. 

Second, it is interesting to probe the holographic model of the nuclear matter phase against the predictions 
of other effective models. In particular, there are claims, most notably from the analysis of baryons lattices 
in Skyrme model \21\ 122) , that chiral symmetry restoring phase transition may take place in the nuclear 
matter phase at zero temperature and finite density 5 . In the skyrmion description this transition is 
also accompanied by the change of symmetry from skyrmion to half-skyrmion [231 124] . Since the Skyrme 
model is naturally embedded in holography, it is important to understand how those results are matched 
by holographic models and ultimately verify the statement about the chiral symmetry restoration (e.g. the 
work of Rho et al. [20 ). If the restoring transition is to take place at all this should definitely happen in the 
nuclear matter phase. 

In this paper we will mostly address the first piece of the above motivation, while restricting to a summary 
of our expectations for the second. We consider two toy-models. In the first toy-model baryons are instantons 
in the limit of zero size, that is point charges in 5D. Since curvature of the fifth dimension binds the charges 
to the bottom of the flavor brane geometry, in the low density phase they form a 3D lattice. However if 
the high enough pressure is applied the 3D lattice becomes unstable and the baryons start occupying higher 
positions in the transverse fifth dimension. More precisely the lattice splits into several sub-lattices separated 
along the transverse dimension. The splitting continues as pressure is further increased. 

In the second toy model we consider a ID chain of baryons, described as finite size instantons in the 
generalized Sakai-Sugimoto model. To render the chain of repelling instantons stable in 3D we turn on 
curvature in all the transverse directions, including the spatial ones. The outcome of the analysis is the 
following. For large instanton spacing the periodic arrangement of instantons takes the form of a one 
dimensional straight chain with an abelian phase difference between neighboring instantons of <j) — 7r. For a 
separation distance, below a critical value given in (|117|) the straight chain becomes an unstable configuration 
and instead the chain takes the form of a zigzag structure (see figure I10[) . The zigzag amplitude e as a 
function of d is drawn in figure [T4l We further find that assuming an abelian phase between adjacent centers, 
the preferred value of the phase is <f> = it for the spacing close to critical. We also consider the case of general 
non-abelian orientation angles and find out that phases with abelian orientations are preferable. We show 
that for small amplitude of the zigzag the neighboring instantons remain antiparallel as in the e = case. 
For larger zigzag amplitude, the relative orientation of the instantons changes from (j) = tt to cj) ~ 117° in 
a first order phase transition. At even larger densities the orientation changes smoothly to an asymptotical 
value 7r/2. For e ^> d the instantons in each of the two layers of the zigzag become closer to each other than 
the original nearest neighbors and prefer to orient themselves in an antiferromagnetic way, <j) — tt. Figure [TBI 
summarizes the phase diagram of the chain up to 2 layers. 

Notice that the approach used in this work is in fact independent of the concrete details of the holographic 
model. Although we start our discussion from Sakai-Sugimoto model and derive our effective action from 
it, the analysis is restricted only to the first non-trivial order in the 1/A expansion. Up to this order all the 
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dependence on a particular model can be absorbed in the definition of the coefficients that determine the 
characteristic scale, such as strong coupling scale A (equivalently the Kaluza-Klein scale Mkk) or 't Hooft 
coupling. Therefore the analysis itself is the same for all holographic models with large A, e.g. the flavored [25] 
Klebanov-Strassler geometry [26] . 

The paper is organized as follows. In section [5] we review the state of the art of the QCD phase diagram 
at finite temperature and chemical potential. Section 12.11 describes realistic QCD with 3 colors and 2- 
3 (massless) flavors. In section 12.21 we summarize possible modifications of the phase diagram when the 
number of colors N c is taken to be large. We also briefly mention there the consequences of large 't Hooft 
coupling A. In section |2~31 we discuss the insight on the phase diagram obtained from the Skyrme model. 

In section |3] we review the story of baryons in holography. For a prototype holographic model of QCD 
we consider a generalization of the one proposed by Sakai and Sugimoto model. Section [3. II is dedicated to a 
review of the relevant facts about the generalization of Sakai-Sugimoto model. The baryons are introduced 
in section 13.21 as instanons of the effective theory on the flavor D-branes. Here we summarize the general 
strategy to be used in our analysis of multi-baryon systems. Section 13.31 contains a brief discussion of other 
holographic models of baryons, which also admit instanton description of baryons. 

Section@]starts with a review of previous results on the phase diagram of holographic QCD. In section R~2l 
we specify the part of the phase diagram that we are going to explore and summarize our expectations. We 
illustrate our expectations through a simple toy model of lattices of point charges, considered in the following 
section [5J In particular in section 15.11 we consider a warm-up example of a ID chain of point charges, while 
section I5~2l generalizes the result to 3D cubic lattice of point charges. 

The effects of finite size are analyzed in section [6l We start in section 16.11 from a review of the BPST 
instanton solution describing an infinite straight chain with a constant orientation twist between neighboring 
pairs of instantons. Later we compute 1/A corrections to the energy of the BPS configuration, which are 
generated by the non-BPS interaction of baryons. Assuming anti-ferromagnetic orientation of instantons in 
section 15721 we show that for large density the straight chain is unstable against splitting along the holographic 
direction (zigzag deformation) in a second order phase transition. The analysis is generalized in section [6~3l to 
the case of arbitrary twist between neighboring instantons. Here we discover another (first order) transition 
at larger densities. 

In section [7] we summarize our results and their interpretation and propose several open problems. We 
supplement the main text by several appendices. Appendix [X] provides the details of calculations of the 
instanton self-interaction energy in the limits of large and small densities. Appendix [B] is dedicated to a 
proof of the fact that the straight chain always prefers anti-ferromagnetic orientation of instantons. Somewhat 
unexpected instability of the zigzag towards formation of the straight chain at high density is demonstrated 
in appendix[C] Finally in appendix|D]we provide the details of the derivation of general twist zigzag instanton 
solution. 

2 Baryons at high density 
2.1 Phase diagram of QCD 

The precise shape of the QCD phase diagram is not firmly established. The current state of the art diagram 
(figure [lj is a collection of lattice results, perturbative and non-perturbative estimates as well as educated 
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Figure 1: Phase diagram of QCD (contemporary view [HSU]). 



guesses at various corners [TJ [2J [3l |4] . In fact not only the boundaries of the many phases of QCD cannot 
be fixed, but even the existence of some of the phases may be uncertain due to a crucial dependence on the 
true values of the parameters, e.g. quark masses, number of flavors etc. 

The low temperature and low chemical potential corner of the diagram is the hadronic phase. In this 
phase quarks are confined into hadrons and chiral symmetry is broken. When the temperature is increased 
a deconfinement transition occurs accompanied by the restoration of chiral symmetry. Let us review here 
both the case of massless and massive quarks. 

In the approximation of massless quarks chiral symmetry is broken spontaneously by the quark conden- 
sate. For high temperature thermal fluctuations destroy the condensate and a symmetry restoring phase 
transition occurs. Lattice simulations in general support the conjecture that at zero chemical potential this 
phase transition is second order for two and first order for three massless flavors. On the other hand at 
finite (large) chemical potential many model calculations suggest that the phase transition to the restored 
symmetry phase becomes first order for both Nf = 2 and Nf = 3. Therefore for the case of two massless 
quarks a change from second order to first order phase transition is expected along the boundary of the 
hadronic and chiral symmetry restored (deconfined) phases. The point where the change occurs is called 
tri-critical (figure [2]). 
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In the case of massive quarks, m u ~ -C m s , the symmetry is broken explicitly. As a result the second 
order transition from a confined to a deconfined phase becomes a crossover, since in the absence of exact 
symmetry there is no critical behavior. Study of QCD at the crossover is extremely complicated, because 
neither the conventional strong coupling nor weak coupling descriptions work in this regime. The phase of 
QCD at the crossover is usually referred as to strongly coupled quark-gluon plasma (sQGP). The crossover 
changes to a first order transition at a (chiral) critical point similar to the critical point of water (figure [1]). 

At small temperature and finite chemical potential physics of QCD is even more interesting. The natural 
energy scale for the interaction of the elementary degrees of freedom (quarks) is the Fermi momentum, which 
is related to the chemical potential as /i 2 = k% + m 2 B . Therefore for large enough chemical potential the 
physics is expected to be weakly coupled. The naive perturbation theory does not work though, because the 
perturbative ground state (free quark Fermi spheres) is unstable due to attraction of the quarks with equal 
and opposite momentum at the Fermi surface. This problem can be overcome in a fashion similar to the 
BCS theory of superconductivity. In the true ground state the quarks form Cooper pairs (condensate) while 
the colored excitations acquire a gap. 

At asymptotically large chemical potential, one can ignore quark masses (with all 3 quark flavors partic- 
ipating). In this case the preferred quark condensate is the one that preserves maximal symmetry: 

<^<&(P)9&(-P)> = ~(^ $ C R (p)<lf R (-p)) = A(p 2 )e ijA e ahA , (1) 

where Greek letters are reserved for spinor indices, a, b and i, j are flavor and color indices respectively. Such 
a condensate breaks both color and flavor symmetries, but leaves the diagonal subgroup of their product 
unbroken. The symmetry breaking proceeds as follows 

SU(3) C x SU(3) L x SU(3) R -> SU(3) CFL ■ 

The color and flavor transformations are locked together, hence the name of the phase: color-flavor-locking 
(CFL). There are in fact two condensates for the left and right fields separately. One of them locks SU(3)l, 
while the other one SU(3)r. So indirectly the chiral symmetry is broken. 

Breaking of the color symmetry results in masses for all gluons. This is analogous to the Meissner effect 
for ordinary superconductors. Therefore the CFL phase is superconducting. Since in the CFL phase all the 
symmetries are broken, while in the QGP they are unbroken, there should be at least one phase transition 
relating these phases. The phase diagram of QCD with three massless quark should thus look as shown in 
figure |2Ja). It is not clear however if there is a phase transition between CFL phase and the nuclear matter 
phase. Since the former phase is deconfined while the latter is confined first order transition can be expected. 
However this phase transition will not restore the chiral symmetry, since it is broken on both sides. The 
separation of deconfincment and chiral symmetry restoration can in principle occur as we discuss in more 
detail below. 

For intermediate values of the chemical potential the difference in the quark masses cannot be ignored. 
A realistic approximation would be to choose = m u = <C m s . First consider infinite m s . In this case a 
condensate of quarks with two out of three colors is formed, e.g. 

<e«3^gi(p)gj(-p)> • (2) 
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The condensate picks a direction in the color space. SU(3) C group breaks down to the SU (2) c and five gluons 
get mass. Gaps are formed for excitations of colors orthogonal to the condensate direction. The potentially 
light degrees of freedom are quarks of the condensate color and gluons of the unbroken subgroup. This phase 
is known as 2SC. In the case of two massless flavors this is the preferable phase at small temperature and 
large chemical potential as shown in figure [3Jb). 

The condensate is invariant under SU(2)l x SU(2)r, which means the 2SC phase has unbroken chiral 
symmetry. Therefore a phase transition must separate this phase from the hadronic phase. This transition 
must be first order, since it involves a competition between the color and chiral condensates. Since there are 
no broken global symmetries the 2SC phase is not a superfluid. 

For finite m s the CFL will be a dominating phase at large chemical potential. However the CFL phase 
requires that condensates that involve the strange quark and a light quark must produce a gap which is larger 
than roughly m^/2/i, which is the difference of the s-quark and light quark Fermi momenta. For not-so-large 
chemical potential the phase is 2SC as shown in figure [T] Since the latter has restored chiral symmetry there 
is a phase transition in between. It cannot be second order, because this will require condensates (us), (ds) 
to be infinitesimal at the phase transition on the CFL side. 

Finally we will discuss one more phase of QCD, which is of most interest for this paper, namely the 
nuclear matter phase. For a chemical potential of order /i ~ lGeV and slightly below the baryons form a 
dilute gas: since it is hard to excite baryons, at zero temperature there is vacuum (a state with no baryons), 
while at finite temperature one expects a small population of baryons in the thermal ensemble. The baryons 
attract each other with a characteristic binding energy Eb, therefore it should cost nothing to produce them 
when \x reaches the value ~ ms — Eb, where tub is the mass of the lightest baryons. The zero temperature 
case is demonstrated by figure J3ja): at small chemical potential the density is zero, while for fi = raj — Ej, 
there is a first order phase transition to the equilibrium density of nuclear matter. If the interaction between 
baryons had been repulsive as in the cases discussed below the transition would have been second order as 
in figure G2h). 




m B -E b m B fj. m B fl 

(a) (b) 

Figure 3: Zero temperature dependence of the baryon density on chemical potential in the case of attractive (a) and 
repulsive (b) equilibrium force between baryons. In the attractive case the transition to non-zero density is a sharp 
first order transition at chemical potential ms — In the repulsive case it is a smooth second order transition. 

In reality this picture is altered by the presence of the Coulomb repulsion between baryons. Because the 
Coulomb force is long range it precludes the existence of the bulk nuclear matter. Therefore the behavior 
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corresponding to figure GUJa) is an idealization, which assumes that the bulk Coulomb force is screened by 
either finite density of electrons or by idealized light s-quark. The real-life QCD is expected to have a first 
order gas- liquid phase transition similar to the phase transition in water. At finite temperature this phase 
transition must end at a (hadronic) critical point, because the Lorentz symmetry will be broken on both 
sides of the phase transition. The end point of the transition should appear at temperatures of the order of 
baryonic binding energy Ej,. 

As follows from the properties of ordinary nuclear matter, it is in a liquid phase. For increased chemical 
potential one may expect the liquid of baryons to turn solid. However it is not clear if such a phase exist at 
all in real life QCD. It may be that the liquid phase goes directly to a dense quark matter phase through a 
deconfining and chiral symmetry restoring phase transition. 

In this part we have reviewed the phase diagram of real-life QCD with N c — 3 and three light flavors of 
quarks. The result is summarized in figure [T] Similar diagrams with more detailed explanation and relevant 
references can be found in the reviews [IJ |3] . 

2.2 QCD at large N c 

Various properties of QCD are often studied qualitatively by taking the limit of large number of colors 
N c |27) . In this limit contribution of quark loops is suppressed and only planar gluon diagrams contribute. 
At zero temperature and chemical potential the large N c QCD is a theory of free mesons and glueballs, 
since there interactions are suppressed by powers of N c , except for interactions with baryons. Masses of 
the mesons and glueballs are of order Aqgd, that is they scale as N®. Baryons are special objects in the 
large N c limit 28; . Because they are made out of N c quarks they are heavier than the mesons and glueballs 
with the mass of order N c Kqcd- Their interaction energy is also non- vanishing and scales as N c . Based on 
these properties of the large N c theory one can speculate on the structure of the T — \i phase diagram (see 
review [5] dedicated to this problem and references therein). One possibility is illustrated by figure [4] Below 
we will briefly explain this diagram. 
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Figure 4: Conjectured phase diagram for QCD at large N c [5] 

At small temperature and chemical potential large N c QCD is a theory of mesons and glueballs. Since the 
baryons are so heavy (~ N c ) their relative abundance is exponentially suppressed. When the temperature 
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is increased a deconfining phase transition is expected. The temperature Td of this transition should be of 
the order Aqcd and the transition is first order. It is implicitly assumed that Nf <C N c . In this regime 
the quark loops are suppressed as Nf/N c and the dynamics is completely captured by gluons, unless the 
chemical potential for quarks (fi q = n/N c ) scales as a power of N c . Thus for /i 9 of order AT? the deconfining 
temperature is expected to be independent from /j, q . 

An interesting question is whether the deconfining phase transition coincides with the chiral symmetry 
restoration or not. Several field theory arguments, e.g. [351 OS, suggest that at fJ, q = Q there is a direct 
relation between confinement and chiral symmetry breaking. However at non-zero chemical potential this 
relation need not hold, and in fact there are various indications to the existence of confined, but chirally 
restored phase [5]. Moreover recent holographic models, which are large N c as well as large 't Hooft coupling 
A models, imply that the two transitions need not coincide even at zero chemical potential [301 . At low 
chemical potential the latter models in fact exhibit phases that are deconfined, but without chiral symmetry. 

For the values of the chemical potential /i above the mass of the lightest baryons (fi q > tub/Nc = m q ) 
and T <Td baryons persist in the nuclear matter phase with a non-zero density. Conversely to the ordinary 
one the large N c nuclear matter must be solid. To see this one can check the ratio of kinetic and potential 
energies in the large N c regime. Potential energy of the baryon-baryon interaction scales like N c . Specifically, 
large N c potential takes a form [3T] 

II ~ N c x A c [r) + N c x A s (r) (Ij ■ I 2 ) (Jr ■ J 2 ) + 

+ iV c x A T (r) (Ii • I 2 ) (3 (n • J x ) (n ■ J 2 ) - (J x ■ J 2 )) + O (1/N C ) , (3) 

where Ac, As and At are the central, spin-spin and tensor potentials independent from N c . Classically this 
potential tries to organize baryons in some kind of a crystal with interatomic separation independent from 
N c . By analogy with the condensed matter systems the kinetic energy is the zero point motion energy in an 
appropriate potential well: 

K ^~±- 1 (4) 

2m B <P N c cPAqcd ' k ' 

where d is the A^-indcpcndcnt diameter of the potential well. Therefore the ratio 

K 1 

n~Nl (5) 

is suppressed in the large N c regime and the cold dense nuclear matter is a crystal. 

The type of the transition between the baryon vacuum and the nuclear matter will depend on whether 
the interaction between baryons is attractive or repulsive. Ordinary nucleons attract each other, but this 
attraction appears to be a result of a very subtle balance of various factors which can easily be violated by 
large Af c effects (e.g. a discussion in (32J [33] .) Thus both of the situations illustrated by figure [3] may be 
realized for large N c . If the baryons still attract the density will behave according to the plot on figure EJa). 
There will be a first order phase transition from the vacuum to the nuclear matter phase at /i = — E^. 
Notice that we leave aside the question of what happens with the Coulomb interaction at large A^ c and its 
effect on nuclear matter. If the baryons repel the transition will be second order, figure G^b). This is the 
situation we will later assume in this paper. 

The crystal phase of cold high density baryons is mostly studied within the framework of Skyrme model. 
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In the next section we will review the main results on Skyrme crystals, while here we continue the discussion 
of the phase diagram. Skyrme crystals exhibit a zero temperature finite density phase transition, which at 
least partially restores chiral symmetry. (This result is also supported by calculations in other qualitative 
models as well [34,.) Cm the other hand it was mentioned before that the deconfining temperature Td should 
be independent of fi when the latter is not too large (quark contribution to the expectation value of the 
Polyakov loop is of order 1/N C , as compared to N® in the deconfined phase). Therefore the chiral phase 
transition and the deconfining phase transition are likely to separate in the dense baryonic phase. 

The authors of [5] refer to the phase in figure [4] below Td and to the right of [i q = m q as "quarkyonic" . 
The reason for this is as follows. Similarly to finite N c QCD at asymptotically large n q a weakly coupled 
description is possible. One can compute the total pressure of finite density QCD perturbatively. The main 
contribution to the pressure is from quarks far from the Fermi surface that scatter with the momentum 
exchange of order fi q . One obtains that the pressure scales as N c and this must be insensitive to non- 
perturbative corrections as one reduces the chemical potential down to [i q ~ Aqc d ■ However for such values 
of \x q the phase should be that of dense baryons also with the pressure of order N c . Effectively this looks like 
the phase is changing smoothly as one changes the chemical potential from (j, q ^> Aqcd down to n q ~ Aqcd 
and both descriptions in terms of quarks and baryons apply. More accurately, for large fi q scattering of 
quarks far from the Fermi surface can be described by perturbation theory. Nevertheless for \x q not too 
large baryons remain the lightest excitations near the Fermi surface: unlike for finite N c QCD, the quarks 
do not screen gluons until the chemical potential fi q scales as a power of N c . The "quarkyonic" phase is a 
phase with quark Fermi sea but with baryonic Fermi surface. When the chemical potential is increased from 
/i q ~ Aqcd the initial baryonic Fermi see smoothly changes to a quark one with the baryonic Fermi surface 
of the width of order Aqcd- 

1 /4 

When the chemical potential /i q is as large as N c Aqcd the contribution of quarks to the pressure 

becomes comparable to that of gluons. Gluons start feeling quarks and the deconfining temperature Td 

l li 

should no longer be independent from [i q . At \i q ~ N c Aqcd the quark contribution completely dominates, 
the gluons are screened and the deconfining transition should end, either at Td = or at some critical point. 

At large N c for /i ^> Aqcd the phase of the quarkyonic matter is perhaps not a color superconductor. 
It was argued in [35] that color singlet chiral density wave (CDW) phase is dominating there. For large but 
finite N c it was proven that the CDW does not dominate until N c ~ 10 3 [36] . 

In holographic models apart from the limit of large number of colors one also has to assume large values 
of the 't Hooft coupling constant A. Since the main topic of this work is a holographic description of the 
nuclear matter phase let us discuss the implication of large A limit on the relevant part of the phase diagram. 
First the baryon mass is proportional to both N c and A. In turn the interaction of baryons is only of order 
A . Therefore baryons are heavy, weakly interacting objects in the large A limit. Specifically the binding 
energy of baryons, if the net force between them is attractive, will be small anyway compared with the 
baryon mass. 

Because of the subleading scaling of the baryon interactions the characteristic scale of baryon physics 
will be proportional to 1/y/X. A characteristic baryon size will be A -1 / 2 and the characteristic density of the 
D-dimensional baryon crystal p oc \ D I 2 . Such densities correspond to the values of the chemical potential 
/i ~ (i c + 0(A ), that is p, — /i c <C /i c - At larger values of the chemical potential, e.g. p — /i c ~ \i c the density 
will be too large for individual baryons to retain their identities. At such densities the nuclear matter will 
turn into some kind of a quark matter. The nuclear matter phase will thus be confined to a narrow window 
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of values of the chemical potential. This is illustrated by figure [5] where we show again both the case of 
attractive and repelling baryons. 




m B fi m B fi 



(a) (b) 

Figure 5: Density as a function of the chemical potential in the large A regime for attracting (a) and repelling (b) 
baryons. The nuclear matter phase is confined to a narrow window of the order Afi ~ ms/A. In the naive diagram 
the transition occurs directly from the no-baryon to a quark phase. 



2.3 High density baryons in the Skyrme model 

In the 1980s the old idea of Skyrme that the baryons can be described as solitons of the chiral model [21] 
received a revived attention |28j . These classical solutions, the skyrmions, satisfy the equations of motion 
derived from the Lagrangian proposed by Skyrme: 

£ = & ttdpUWrf + ^tr %UU\d v UU^ 2 , (6) 

where U (x) = a + fr(x) ■ f is the SU (2)-vahied pion field and f n is the pion decay constant. The second term 
in the Skyrme Lagrangian does not follow from chiral Lagrangian, but was added by Skyrme to stabilize the 
solutions. The parameter e is typically fixed such that the mass of the single skyrmion matches the masses 
of N and A nucleons. Skyrme model is effectively a model of large iV c baryons. 

Crystals of skyrmions were extensively studied in the 80's starting from the work of Klebanov [25]. To 
find a skyrmion solution corresponding to a lattice one imposes appropriate periodicity conditions on a single 
skyrmion in the unit cell of the lattice, consistent with the symmetries of the given lattice. The solution for 
various types of cubic lattices is typically obtained numerically. Analytical results are possible only within 
an approximation scheme, e.g. replacing the unit cube by a 3-sphere |37j . or employing the Atiyah-Manton 
ansatz [35]. (See also [39] for other approximation techniques.) 

It was noticed that at higher densities skyrmion solutions exhibit an extended symmetry as compared 
with the symmetries allowed by the original lattice |40[ 141] . This is most obvious in the S 1 -approximation 
of the cubic cell. When the size of the spherically symmetric skyrmion is small compared with the volume 
of S 3 it has only an SO (3) symmetry. However when its size becomes comparable to that of the sphere 
the symmetry extends to 50(4). In the case of simple cubic lattice the additional symmetry corresponds 
to an exchange of the vertices of the elementary cells of the lattice with their centers. Instead of the cubic 
lattice it becomes more natural to consider unit cells of half the original volume consistent with the extended 
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symmetry (Wigner-Seitz cells of the bcc lattice). The new symmetry imposes strong constraints on the 
skyrmion solution. As was demonstrated by Goldhaber and Manton 23J each elementary cell contains a 
half of the skyrmion charge and the value of the skyrmion field is different for the vertices of the original 
cube and its center (in particular it differs by the sign of a). Thus instead of the original cubic lattice of 
skyrmions one obtains a bcc lattice of half skyrmions with different kind of skyrmions in the sites and in the 
centers of the cubic cell. It was also shown that, when averaged over the original unit cell, one obtains 



where brackets mean average over space, since the Skyrme model is classical. This serves as an indication 
of the chiral symmetry restoration at high density and zero temperature consistent with the phase diagram 
in figure 2] 

It is not so easy to see chiral symmetry restoration in the case of cubic lattice, but it becomes obvious 
if one resorts to the spherical cell approximation [42]. When the skyrmion on S 3 is small it necessarily 
points in some direction in the group space. The solution is roughly a map R 3 — > S 3 as in the case of the 
flat-space skyrmion, with U — > 1 at inifinity. However, when the size of the skyrmion is of order of the 
S 3 size the solution becomes the identity map S 3 — > S 3 of the spatial 3-sphere onto the group space. The 
integral of U over the sphere vanishes in this case and specifically (O holds. Further studies of excitations 
over the half-skyrmion lattice have shown the doubling of the spectrum, consistent with chiral symmetry 
restoration [4"5] . 

It was shown by Kugler and Shtrikman [24] that the preferred configuration of skyrmions at small density 
is the fee lattice, for which each skyrmion has 12 nearest neighbors. At higher densities this lattice becomes 
a simple cubic lattice of half skyrmions. The sc lattice of half-skyrmions can attain the smallest possible 
energy per skyrmion out of all configuration of skyrmions. This happens at the density corresponding to the 
lattice spacing d = 4.7 in units of l/e/ ff . This gives the energy E — 1.038 in units of the topological lower 
bound on the skyrmion energy (to be compared with E = 1.23 for a single skyrmion). As in the case of the 
Goldhaber-Manton transition the transition in the fee skyrmion lattice also points to the chiral symmetry 
restoration. 

It turns out that skyrmions can naturally be related to the baryons in holographic models. As will be 
discussed below baryons in holography are described by instanton solutions of the gauge fields on flavor 
branes. Remarkably the holographic description of baryons by instantons was anticipated by Atiyah and 
Manton [38], who suggested the following approximation for skyrmions. The skyrmion field U in 3D can be 
approximated by the holonomy of a 4D Euclidean instanton field. More precisely 



where the holonomy is evaluated for the additional 4th component of the gauge field along the additional 
direction. The instanton number in this case is mapped to the skyrmion number. It was later realized that 
such an approximation works extremely well. In particular, for a single skyrmion it predicts the energy 
E = 1.24, while for the minimum energy of the Kugler-Shtrikman lattice of half-skyrmions the energy 
predicted by Atiyah-Manton ansatz is E = 1.04. 

For a while the coincidence of the results from skyrmion and instanton pictures remained somewhat 



(a) = (tt) = 



(7) 




(8) 
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mysterious. It was in fact holography which provided a resolution for this mystery [44] . In some holographic 
models baryons are described by instanton solutions of the non-abelian gauge fields that live on D-branes, 
e.g. Sakai-Sugimoto model [?:. For static baryons the instanton solution lives in 3 spatial and the holographic 
dimension "z" . The pion-skyrmion held U naturally appears in the holographic setting as the holonomy of 
the component of the non-abelian field along the holographic direction, cf. ([5]), 



The Skyrme Lagrangian ([6]) as seen by holography will be modified. An infinite number of terms coupling 
the pion field to the infinite tower of vector mesons will be added to it. This "improved" Skyrme model was 
derived by Nawa, Suganuma and Kojo in [45] and further used to study various properties of baryons. 

So why does instanton approximation describe skyrmion solutions so well? Holography provides for- 
mula §§§ for the skyrmion of the improved model. Unlike the approximation of Atiyah and Manton the 
instanon of equation ([9]) lives in a curved space. The curved space is in fact responsible for the tower of 
vector mesons that couple to the skyrmion. In the limit of large 't Hooft coupling A the instanton size is 
much smaller than the curvature scale and one can think of it as of just a flat space instanton. Thus the 
formula of Atiyah and Manton gives an approximation to the improved Skyrme model and further to its 
truncation, the original Skyrme model. The fact that the approximation is so good can be attributed to the 
large value of 't Hooft coupling constant on the holographic side. 

With the straightforward relation between skyrmions and holographic baryons a natural question to 
address is how the results on the skyrmion lattices can be translated to a holographic setup. It is also 
interesting if the chiral symmetry restoration at zero temperature and high density can be proved using 
holography. A qualitative proposal on mapping the skyrmion results to holography was made in the paper 
by Rho, Sin and Zahed |20j . Later by studying toy models with ID baryon lattices we will find some support 
to the ideas of that proposal. It still remains an open problem though to give an explicit description of 
skyrmion physics in the 3D case. Now let us move to a discussion of baryons in holographic models. 

3 Baryons in the Sakai-Sugimoto model 
3.1 The model 

The holographic model proposed by Sakai and Sugimoto [7] describes large N c , large A QCD in the quenched 
approximation Nf <C N c . Geometrically Sakai-Sugimoto model corresponds to an extension of the Witten's 
model [5] through introducing Nf probe D8 and anti-D8 branes. In the Witten's model one considers N c 
type IIA D4-branes in the near-horizon limit. Four of the dimensions transverse to the D4 are compactified 
on the 4-sphere S 4 . One dimension along the D4 is compactified on the circle S 1 with antiperiodic boundary 
conditions for fermions. The latter procedure breaks supersymmetry and introduces a (Kaluza-Klein) scale 
Mkk into the theory, which is related to the radius R of S 1 . Correspondingly the geometry breaks into a 
warped product 



For such a geometry to be smooth the radial direction R + should end at a finite distance from the origin, 
call it ua following the original papers. The value of ma is controlled by the Mkk scale or R. The R + x S 1 




(9) 



AdS 6 xS 4 4 M 1 ' 3 x R+ x S 1 x S 4 . 
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part of the product thus looks like a cigar as shown in figure [51 

The geometry of the D4-brane is described by the type IIA supergravity solution for the metric, the 
4-form and the dilaton 
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where the curvature radius of AcISq and warping f(u) is given by 

Rli = Kg s NJ 3 sl f(u) 



i 



(t)' 



(10) 



(11) 



V4 is the volume of the unit 4-sphere, l s — \fa! is the string length and g s is the string coupling. The radial 
coordinate u runs from ua to infinity. In terms of the above metric the relation between u\ and the circle 
radius R reads 

*>D4 



2nR = T 



at: rm xl/2 



(12) 



The theory on the world- volume of the D4-branes is a 5D Yang-Mills theory compactified on 5 1 . An- 
tiperiodic boundary conditions break supersymmetry: the fermions receive mass of order Mkk- In the low 
energy limit the theory is a pure glue SU{N C ) Yang-Mills with the strong coupling scale A set by Mkk- 

From the DBI action for the D4-branes one reads the coupling of the 5D Yang-Mills theory [10] 

gl = (2tt) 2 9s I s . (13) 

Below the compactification scale Mkk the theory is effectively 4-dimensional with the coupling 

g 2 =g 2 5 /(27rR). (14) 

Note that the theory has a running dilaton so the effective coupling g\ oc diverges in the UV limit, u — > 00. 
Although the theory constructed so far is not UV-complete, one should not worry about this when the IR 
properties are studied. The 't Hooft coupling A is introduced as 



A = glN c . 



(15) 



To introduce flavor degrees of freedom into the game one immerses Nf D8 and anti-D8 flavor branes into 
the geometry of D4 [7]. In such a scenario strings that connect D4 and D8 (anti-D8) carry quark (anti-quark) 
degrees of freedom. In general solutions for the geometry of D4-D8 systems are not known: one typically 
resorts to the probe approximation, in which the backreaction of the flavor branes on the geometry of D4 is 
ignored. The latter regime is achieved when Nf <C N c . 

In the setup of the Sakai-Sugimoto model the D8-branes span the Minkowski directions and wrap the 
S in the transverse directions. The remaining dimension spans a ID curve on the cigar R + x S . At zero 
temperature stable embeddings correspond to the D8-anti-D8 connecting smoothly in a U-shape configuration 
(figure [S]). This is a geometrical realization of the chiral symmetry breaking. By default there are two gauge 
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Figure 6: The figure on the right is the generalized non-antipodal configuration. The figure on the right 
describes the limiting antipodal case L — 7r_R, where the branes connect at uq — u\. 

groups U(Nf)z, and U(Nf)R which live on the D8 and anti-D8. However when the D8 and anti-D8 reconnect 
only the diagonal subgroup of the product survives: 

U(N f ) L x U(N f ) R -> U(N f ). 

The U-shape embcddings are parameterized by the asymptotic separation of the D8-anti-D8 branes on 
the circle - L, or equivalently by uq - the value of the it-coordinate of the embedding closest to the origin 
(figure^. For uo = u\ the configuration is called antipodal since the D8 and anti-D8 remain at the opposite 
points on the circle L = ttR for all values of it. Note that in general there is no parameter in real QCD 
analogous to this new scale L (ito). 

The background configuration of the Sakai-Sugimoto model is confining. The low-energy degrees of 
freedom are glueballs (closed strings), mesons (strings ending on both D8 and anti-D8) and baryons (N c 
strings ending at a point in the bulk). If one puts the Sakai-Sugimoto model at finite temperature - there is 
a deconfining phase transition |30) . At finite temperature the Euclidean time direction is also compactificd 
and the background corresponds to a black hole. When the temperature is increased to the critical value 
T = 1/2ttR the geometry changes: the radial direction R + and the Euclidean time now have the topology 
of a cigar, while R + x S 1 becomes a cylinder. Such a swap corresponds to a first order phase transition. 
There is a number of ways to see that this is the deconfining phase transition, e.g. one can look at the 
quark-antiquark potential or scaling of the free energy with the number of colors N c in the two phases. 

Depending on the value of L one can have a chiral symmetry restoration either together with deconfinc- 
ment or at some large value of the temperature: two transitions do not necessarily coincide. Geometrically 
chiral symmetry restored phase corresponds to the D8 and anti-D8 separating from each other. In the large 
temperature phase the D8 branes live on the cylinder rather than on the cigar. Since there is a horizon 
(minimal value of the it-coordinate on the thermal cigar) , there is a natural place for the D8 to end. There 
exist two classes of embeddings: the U-shape and parallel D8-anti-D8-branes ending under the horizon. At 
high enough temperature the second class becomes favorable. 

The low-energy dynamics of the flavored degrees of freedom is described by the effective action comprised 
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of the Dirac-Born-Infcld (DBI) action and the Chern- Simons (CS) term: 



S = Sdbi + s cs ■ (16) 

The DBI term reads! 

Sdbi = T 8 J d 9 xe^str - det{g mn + 2™' J" mn )) , (17) 

D8+D8 

where T p — (2tt)^ p 1^ p ^ 1 is the D8-brane's tension, g mn is the nine-dimensional induced metric on the branes. 
Fmn is the U(Nf) gauge field strength on the world- volume of the D8-branes, 

T = dA + Af\A, (18) 

and "str" is the symmetrized trace over the flavor indices. The CS term comes from the coupling of the 
D8-branes to the bulk RR field, 

Scs=T 8 J C 3 Atre 2 ™'^, where F A = dC 3 . (19) 

D8+D8 

In terms of the effective action low energy degrees of freedom correspond to fluctuations of the embedding 
profile and of the gauge fields. In the weak field approximation (weak compared to the scale set by A or 
equivalently curvature) one can restrict to the leading orders in powers of T mn in the action: 

Sdbi = So [g] + S YM + (T 4 ) , (20) 

where Sq is the action with the gauge fields set to zero, while <SVm corresponds to the first non-trivial order 
in gauge fields. The low-energy modes are independent from the coordinates on S 4 : one can integrate out 
this directions and work with an effective 5D theory. Since the gauge fields along the S 4 directions will not 
play any role in our latter discussion and will be switched off, the Sym becomes a 5D U(Nt) Yang-Mills 
theory with a coupling depending on the holographic coordinate. In this respect a more natural way to 
parameterize the holographic direction is through a new coordinate, for which the action simply takes the 
form 

where the indices are raised with rf' v . The 5D Yang-Mills coupling constant is given by the equation 



Mkk—, where — = J . (22) 



9ym( z ) "a dz V u 5 R 3 L 



D4 



The values of z run from — oo to +oo, with positive (negative) branch corresponding to the D8 (anti-D8) 







2 We reserve the lowercase Latin indices for the world-volume of the D8-branes, Greek indices will parameterize coordinates 
in the Minkowski space and holographic direction, while the uppercase Latin indices will label four spatial indices in the latter 
5D space-time. 

3 Note that the coordinate z here is not the same as in case of the original work of Sakai and Sugimoto [7] . It rather coincides 
with w introduced in 1321 . In terms of z coordinate all components of the 5D metric are the same. 
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The coefficient in front of the action is often denoted as K. Using (fT2"j) and (|15p it can be written in terms 
of the 't Hooft coupling and N c 



K = n77~~3 ■ (23) 



AiV c 

2167T 3 

For our analysis below there is no need to know the explicit form of u{z). In fact all we need is the expansion 
of u for small z 

- = c(l + ^^M* KK z 2 + 0(M* KZ *)), C=^- (24) 

Finally we are interested in the case of two flavors Nf = 2. It is convenient to decompose the U(2) held 
A into SU(2) and U(l) parts 

A = A su(2) + i A u(1) . (25) 
The effective action of the Nt = 2 Sakai-Sugimoto reads 

S[A,A] =S YM + S CS , (26) 

where Yang-Mills action (|2Tj) takes the following form in terms of the U(l) and SU(2) fields: 

S Y u = -* J d** d, ^_ (tr Fl + \ J*) , (27) 

The CS term (fT9")) integrated over the S 4 reduces to 

s - = 2&A(-^-^-^ 5 )' (28 » 



For iVf = 2 the only terms that remain are 



S CS = -^ ( AMvF 2 + /iAF 2 . (29) 



16tt 2 7 96tt 2 

The above derivation ignores backreaction of the gauge fields on the geometry of embedding. In fact 
non-trivial values of the gauge fields on the D8-branes may change the embedding, e.g. [IT] ■ In the weak 
field approximation the backreaction can be taken into account by adding scalar fields describing transverse 
fluctuations of the brane profile. The effective action for these scalar fields was derived in 32 . In particular 
it was shown that the leading order contribution comes from an isosinglet scalar field $. Introduction of the 
scalar merely results in a renormalization of the effects induced by the abelian field A. But before discussing 
these effects in detail let us remind how this effective action is linked to baryons. 

3.2 Baryons in the Sakai-Sugimoto model 

To introduce baryons in holography one has to have an object, where N c quark strings can end - the baryon 
vertex (b.v.) [TJ]. One way to construct it is to consider Dp-branes wrapped on compact p-cycles with N c 
units of flux. In the above setup this works as follows. There are iV c units of the i^-flux through S* 4 , 



1 



(2^ 



Fi = N c . (30) 
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Introduce a D4-brane wrapped on S* 4 . The D4 is coupled to the RR field. As a result it picks N c units of 
charge from the U(l) gauge field B on its world- volume: 



Due to the charge conservation the net charge in a closed space must vanish. To compensate for N c units of 
charge the same number of strings must end on the wrapped D4, since each of them carries precisely one unit 
of charge. This is an appropriate extension of a b.v. construction proposed by Witten for AdSs x S 5 [12] 
and its generalization to confining backgrounds [4"5] . 

In the non-supersymmetric background of Sakai and Sugimoto the strings pull the b.v. towards the D8- 
branes. For all values of uq (L, figure the minimum energy configuration corresponds to the b.v. sitting 
precisely at uq [471 . In this case the b.v. "dissolves" in the D8-branes and can be described by an instanton 
- SU(Nf) flavor gauge field configuration in their world- volume [48] . 

It is clear from the Chern-Simons action (f2T)|) that the U (1) field A is sourced by the topological density of 
instantons. Indeed since we deal with the vector-like t/(l) of the product U(Nf)L X U(Nj)r spontaneously 
broken by the U-shape configuration, the bulk field A sources the baryon number current in the gauge theory 
on the boundary. According to the standard lore the subleading asymptotic in the z — > oo expansion of A 
corresponds to the vev of the dual operator, baryon charge density (/o_b), and the leading one is the source, 
i.e. the baryon chemical potential fis- Strings ending on the flavor branes (quarks) carry a unit charge under 
this U(l). For a static iV-instanton configuration of the SU(2) fields one obtains from the first term in (|29l) : 



Since the baryon number is 1/N C times the U(l) charge, for the particular normalization of the U(l) generator 
in (|25|) it is a half of the instanton number N. 

We will be interested in static baryon configurations, that is the instanton-like solutions for the non- 
abelian fields A^ on the flavor branes. To find equilibrium configurations we must compute the free energy 
on a generic instanton solution and minimize it with respect to parameters, such as position, size and 
orientation. It turns out that the natural size for the instantons is set by the scale 1/MkkvX. Although full 
instanton solutions for the curved space action (|26[) are not known, the above scaling of the instanton radius 
allows to consistently expand the action in powers of 1/A as first analyzed by Hata et al in the work [13] . 
The leading O(A) term in the full action ([26]) is simply the fiat-space 5D Yang-Mills action, <7ym(z) = const. 
The Chern-Simons term is explicitly O(A ), while (AIkkz) 271 curvature corrections to 1/sym wm be of order 
0(A 1_n ). Indeed the n— moment of the instanton density should scale as a™, where a is the radius. 

Thus the only O(A) term in the action (f26| is the SU(Nf) Yang-Mills term with constant coupling. If the 
0(1) terms in the action are ignored, static baryon is a BPS (BPST) instantonQ The instanton configuration 
is given by four Am components of the 5D field, while Aq and A^ vanish. Parameters of the BPS instanton 
solution are free moduli. In other words the energy of a configuration does not depend on the values of the 
instanton parameters at this order. At the O(l) order the interaction terms turn on and (partially) fix the 
moduli. In particular the instanton density sources a notrivial solution for A$, while Aq and Am can still 

4 It can in fact be shown that the zero-size (anti-) self-dual gauge fields localized at z = minimize the full DBI action 




(31) 




(32) 



[13] [25] [32] . 
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be kept trivial. Notice that there is no need to compute corrections to the instanton solution Am at the 
next-to-leading order. Indeed corrections to the equations of motion due to 8 Am will be proportional to 
the variation of the action on the solution itself, which vanishes identically. Such corrections will only affect 
0(A _1 ) terms in the action. Therefore up to the O(l) order the Am fields can be considered as self-dual. 
The O(l) equation for A reads 

nA = -—£±—i(x), (33) 

where the Laplacian is taken in the 4D space of x M and we have introduced the topological density 

I(x) = ^e MNPQ trF MN F p Q. (34) 

As can be inferred from [32] the equation for the isosinglet scalar field 3>, responsible for transverse 
fluctuations of the D8-branes, satisfies a similar equation to (|53")) , except that the scalar couples to tiF MN 
instead of the topological density. However according to the above remark at the O(l) order in 1/A the 
instanton solution is approximately self-dual. Therefore at this order $ satisfies the same equation as Ao up 
to a numerical coefficient. Using (|33[) one can derive the following expression for the baryon free energy at 
the O(l) order: 

^[moduli] = k(M kk Jd 3 xdz (l + M 2 KK z 2 + O (M^z 4 )) \ ^F 2 MN + 

' 1 d 3 xdz A x I{x) + 0(1/ X), (35) 

where the integrals are evaluated on the O(A) order flat space solutions, which define Aq through equa- 
tion ([33]) . Coefficient C takes into account contribution of the scalar field. Since the scalar satisfies equation 
similar to Q33[) it merely rcnormalizes the instanton self interaction term 

1 - C 3 

£7 = 1 . (36) 

For the antipodal configuration of the D8-branes (( = 1) scalar fields have no effect on the profile and the 
scalar field decouples. Conversely the minimum value C can take is 8/9, when C, —¥ oo. In other words for any 
( the coefficient C is positive. This means that in the Sakai-Sugimoto model the short distance interaction 
between baryons is always repulsive, never mind the value of ( [32]. At large distances the sign of the 
interaction is determined by the exponentially decaying Yukawa tails. Since the scalar meson mediating the 
attraction is always heavier than the vector in the Sakai-Sugimoto model, the net baryon force is repulsive 
also at large distances!! 

Expression ()35j) can be recast in a nicer form through a rescaling of parameters: 

8C3_5 M|k->M 2 , ^£_^AVC, io^^- (37) 



Upon this rescaling coefficient \/C becomes an overall multiplicative factor in the expression for energy and 



5 One may also recall that massless pions are also present in the Sakai-Sugimoto model. Therefore they should dominate the 
long-distance force. However in the bulk nuclear matter their contribution is expected to largely average to zero and isoscalar 
exchange is considered dominating. We thank the referee for bringing this issue to our attcnsion. 
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will be dropped from here on. Equation Q35p then takes form 



where 



E [moduli] 



5ym 



-Ena + Eq — / c 
N C XM (l + M 2 z 2 




I 



d 3 xdz A x I{x) + O^ 1 ) , 



(38) 



(39) 



and 




(40) 



Equation (|38l) will be the principle formula for our study of baryon systems below. We will refer to the first 
integral in (|38p as to "non-abelian" term, while the second term will be called "Coulomb" or "abelian" . After 
rescaling (|37[) the new expression for the free energy is quite universal. The rescaling hides model-dependent 
coefficients in the definition of A and M. Therefore we expect (|3"8")) to be correct beyond the Sakai-Sugimoto 
model. Note also that this expression assume net repulsion between baryons (instantons). The adjustment 
for attracting baryons can be achieved e.g. by changing the sign in ([3D]). 

Notice that the expression for the free energy depends on moduli of the instanton solution. This depen- 
dence shows up at O(l) order due to the curvature correction and Coulomb term. These terms stabilize 
instanton moduli. For example in a single instanton case the energy will depend on the z-locus of the in- 
stanton and its size. Curvature terms will prevent instanton from being large and pin it to the bottom of 
the brane configuration z = 0. On the other hand Coulomb self-interaction dislikes small charges and will 
stabilize instanton size at a finite value. 

Apart from the two above parameters, the full moduli space of the flat-space static instanton also contains 
the coordinates in the remaining three spatial dimensions and an SU (2) orientation up to an overall sign. 
For a single instanton these are the true moduli even when the NLO corrections to the energy are taken into 
account. Upon quantization the 5*0(4) rotational symmetry of the orientation moduli space gives rise to the 
spin and isospin quantum numbers, which are locked together I — J. Fermionic baryons have half-integer 
spin (isospin). Spectrum of light baryons was studied in [T3] and [4"T] . 

In this work we will consider classical baryons with a fixed isospin orientation, which are neither fermions 
nor bosons. We will ignore the effects of quantum spin and isospin interactions. Let us use equation (|38p 
to calculate the mass and the size of a single instanton. Evaluating (f3"8")) on the single instanton solution 
localized at z = one finds 



where a is the size of the instanton. Minimizing with respect to a one finds the equilibrium radius of a 
standalone instanton: 




(41) 



(2/5) V 4 
Ma/A 




(42) 



The mass is given by the non-abelian energy. Substituting a for a gives 




(43) 
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As expected the mass of the baryon is proportional to N C A: baryons are heavy in the large N c limit. In 
holography they are even heavier because of large A. 

Now it is clear that the curvature correction to the flat space Yang-Mills action as well as the Chern- 
Simons term give O(A ) contribution. Indeed the (zM) 2n corrections should be proportional to Aoq™ oc A 1- ™, 
because M effectively is the only scale in the problem. The instanton topological density sources the abelian 
field Aq in (|2^|) . From the equation of motion for A (|4T))) it is obvious that the abelian field is subleading to 
the instanton solution. Therefore apart from the O(A) Yang-Mills term in (|2"T|) we must only keep the flat 
space Maxwell action, which is O(l), and drop the second term in the Chern-Simons action (|29|). which is 
only 0(A _1 ). 

An important comment is in order. The size of the real-life baryons is approximately as ~ 4i?Yukawa, 
where R Yukawa must correspond to MZ K . In the large A limit the size acquires additional 1/vA suppres- 
sion (|42p. This signifies that in general it is not sufficient to just rely on the low energy DBI action (|17l) . 
Stringy corrections including derivative terms may be equally important. Nevertheless we hope that the 
higher order corrections will not affect any of the conclusions below, but rather renormalize the exact values 
of physical quantities. 

To summarize the strategy for studying static classical baryons is as follows. First one has to find a 
flat space instanton solution relevant for the problem. In general this is done with the help of ADHM 
construction |49j . A particular formulation of the ADHM construction convenient for baryon study will 
be discussed in section [6] Next one has to solve equation (|40[) for Aq and evaluate the energy (|38|) as a 
function of instanton moduli. The energy then needs to be minimized with respect to moduli to find stable 
configurations. 

In the applications below we will be considering toy models of one dimensional baryon configurations. 
We will typically choose to label the direction along the configuration as X4, while the transverse directions 
will be labeled x\, xi and X3 = z. To stabilize ID configurations it will be useful to introduce an ad hoc 
curvature in all dimensions transverse to X4. In particular we will consider two situations: an isotropic 
transverse curvature, for which the curvature is the same in all transverse directions, 



8tt 
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= N c XM(l + M 2 (x 2 1 +x 2 +x 2 ))- (44) 



and the anisotropic one 



8^ 



= N C XM (l + M 2 x\ + M' 2 (xf + xl)) . (45) 
9ym ^ ' 

Notice that equations (14"2"T) and ([4"5]) will be modified for the above choices. In particular, the equilibrium 

instanton size will be 

(2/!5) 1/4 

isotropic a = ■=— , (46) 

MyA 

(1/5) 1/4 , \ 

anisotropic cto = = , (47) 

y/XMM' 

where in the latter case we have also assumed high anisotropy M' ^> M . 
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3.3 Baryons as instantons in other holographic models 

As was discussed in the previous section holographic nuclear physics based on the generalized Sakai Sugimoto 
model suffers from several drawbacks. The most important ones are: (i) The size of the baryon scales as 
1 / y/X (|42p and hence in the large A limit stringy corrections cannot be faithfully neglected, (ii) Due to the 
fact that the lightest scalar particle that couples to the baryon is heavier than the lightest vector meson the 
interaction between two baryons is dominated in the far zone by repulsion. Thus no nuclei bound states of 
nucleons will be stable in this model. 

This situation naturally raises the issue of constructing other holographic laboratories for nucleons and 
nuclear interactions. Such holographic models should be dual to a boundary field theory that admits confine- 
ment and chiral symmetry which is spontaneously broken. An additional constraint on the construction is 
that the corresponding baryonic vertices will be immersed in the flavor branes such that the baryons can be 
described as instantons of the flavor gauge symmetry. Geometrical realization of the spontaneous breaking of 
flavor chiral symmetry can be achieved, in a similar manner to the Sakai Sugimoto model, by incorporating 
a U shape structure for a stack of Nj branes and anti-branes (see figure [6]) . Indeed such a construction was 
implemented in other models and in particular in the non-critical six dimensional holographic model [50] as 
well as in the model based on introducing D7 anti-D7-branes into the KS model which is described in [25] 
and |5lJR Let us now describe some features of the models and of the corresponding instantons as baryons 
in those models. 



3.3.1 The six dimensional non-critical model 

The non-critical holographic model of [47] is based on six dimensional background associated with large N c 
D4-branes. By compactifying one world-volume coordinate of the D4-branes, the background is rendered 
into a confining one. One then incorporates Nf flavor D4 and anti-D4-branes which are perpendicular to 
the "colored" D4 branes. It is well known that this model, similar to other non-critical models suffers from 
the problem that the curvature of the background cannot be made small by taking the large N c , and it is of 
order unity. 

In the absence of compactificd extra dimensions (in addition to the circle mentioned above) , the baryonic 
vertex in this model is not a wrapped brane but rather simply a DO-brane. One can argue that similarly 
to ([32]) the action on the worldline of the DO brane has a CS term that forces the baryonic vertex to have N c 
strings attached to it and connecting it to the flavor branes. Due to the fact that the model is at A ^ 1, the 
problem of the stringy scale baryon radius does not occur for this model as was shown in [47) . Non-critical 
models also correctly predict certain mo del- independent baryon properties of large N c QCD [53] Q 

The disadvantage of having 't Hooft parameter of order one is in fact advantageous also for the meson 
spectrum. Generically in critical holography the masses of the scalar and vector mesons scale like the scale 
of the system [M ~ 1 /R where R is the radius of the compact direction) and the masses of the higher spin 
mesons, which are described by semi-classical stringy configurations scale as VT ~ M\f\. This implies that 
there is a huge gap between the low and high spin mesons. In the non-critical holographic model on the other 
hand since A ~ 1 there is no such a gap. This feature may be an asset also in terms of the baryon interaction. 

6 We do not attempt to review here other holography inspired models, where baryon physics is intensively studied, e.g. |45II52] , 
In those bottom-up models there is no A, it is rather explicitly or implicitly set to something of order one. Thus the problem 
of large A is replaced by a lack of microscopic or top-down justification of the models. 

7 In the earlier versions of this paper, citing [531 . we mentioned some issues with the derivation of the model-independent 
properties in critical string models. This issue was later resolveld in [54) . We thank the referee for communicating this result. 
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If the mesons that couple to the baryon are stringy mesons that reside on a Regge trajectory, than indeed 
the lightest scalar meson is lighter than the lightest vector meson which will guarantee the domination of 
the attraction over the repulsion in the far zone. 



3.3.2 The DKS model 

This critical model is based on incorporating D7 and anti-D7 flavor branes in the Klebanov-Strassler back- 
ground [33]. It was proposed by Dymarsky, Kuperstein and Sonnenschein (DKS) in [33] (see also an earlier 
realization [55] for the Klebanov-Witten background [55]). The D7-branes admit a U-shape profile whereby 
the spontaneous breaking of the chiral flavor symmetry is geometrically realized. 
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Figure 7: Baryon in the DKS model. For highly non-antipodal embedding Uo 3> Ua baryon is represented by 
the baryon vertex and N c strings stretching towards the D7-branes (a) For nearly-antipodal Uo — Ua configuration 
baryon vertex "dissolves" in the D7-branes and can be described as a world-volume instanton configuration of flavor 
fields (b). 



A baryon in this model consists of a baryonic vertex in the form of a D3-brane wrapping the S 3 of the 
base of the conifold with N c strings connecting it to the flavor branes. Depending on the way the U shape 
flavor branes are embedded in the background, there are two possibilities for the location of the the baryonic 
vertex |51j . When the tip of the U-shape is very far from the tip of the conifold the baryonic vertex will 
approach the tip of the cone, namely, will be separated from the flavor branes, whereas when the tip of the 
U-shape is close to the tip of the conifold the baryonic vertex dissolves in the flavor branes (see figure [7]). In 
this latter case one can describe baryons in terms of flavor instantons in a similar manner that was described 
above for the generalized Sakai-Sugimoto model. 

The leading order interaction between the baryons in the far zone where the separation distance is much 
larger than the inverse of the basic mass scale (the analog of Mkk) is determined via a competition between 
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the attraction and the repulsion due to the exchange of isoscalar scalar mesons ++ and isoscalar vector 
mesons 1 respectively. The unique property of the DKS model is that the lightest ++ particle is in 
fact a pseudo-Goldstone boson associated with the spontaneous breaking of scale invariance. In a certain 
range of parameters this meson is parametrically lighter than any other massive state. As a result unlike the 
Sakai-Sugimoto model, in the DKS model the baryons attract each other in the far zone. At short distances 
the potential admits a repulsive hardcore. 

Thus we see that the two major problems of the Sakai Sugimoto model, namely the size of the baryon and 
its interactions, can be bypassed in other models. The former in the non-critical model and the latter in the 
DKS model. It is still to be seen if the non-critical model admits the desired property of the domination of 
the attraction over the repulsion. Alternatively, a new consistent holographic model where the two problems 
are been solved simultaneously is obviously very welcome. 

In the following sections we discuss the structure of the holographic nuclear matter based on determining 
the lowest energy configurations of the multi instanton solutions. It is important to emphasize that, though 
the analysis was based on the Sakai Sugimoto as a prototype model, it could be adopted also for the cases 
of the non-critical and the DKS model and in fact to any holographic model in which flavor instantons play 
the role of baryons. 

4 Baryonic popcorn 

In this section we turn on finite density for holographic baryons. Due to repulsive interaction between baryons 
it is natural to expect that at large density the condition 2 = (location at the tip of the U shape) for the 
baryons may no longer minimize the energy. In other words one can expect that densely packed baryons 
will overcome the gravity force, pinning them down to the tip, and expand into the holographic dimension. 
Such effect was observed in |18) . where a uniform distribution of instantons was assumed. Formation of a 
macroscopic layer of baryons in the holographic dimension was identified with a developing Fermi sea of 
quarks. Uniform density of baryons may indeed be a good approximation at large enough densities, when 
matter is expected to be quarkyonic. Instead we consider discrete instanton solutions, which should be a 
better description for the nuclear matter phase in the vicinity of phase transition. It is in this regime of 
relatively small densities we expect to find the onset of the quarkyonic phase. 

4.1 Multi-baryon systems 

Sakai-Sugimoto model at finite baryon density was originally considered in |57j . For the analysis of properties 
of dense hadronic matter Kim, Sin and Zahed employed the fact that the asymptotic value of the time 
component Aq of the U(l) gauge field on the world-volume of the D8-branes corresponds to the baryon 
chemical potential in the boundary gauge theory. Phase diagram of large N c QCD was later studied in JTSJ 
[TCI [T71 \TE\ using Sakai-Sugimoto model as a prototype. For other holographic models with finite baryon 
density see e.g. [19]. 

Phase diagram of Sakai-Sugimoto model at finite temperature and density was summarized by Bergman 
et al. in [T7]. In accordance with the review in section [2] phase diagram of the holographic analog of QCD 
has the following features. At low temperature and chemical potential the phase is confined with broken 
chiral symmetry. For the values of chemical potential up until the mass of lightest baryons m b the thermal 



2G 




log77l B '~~ Region of interest 



log ]i 



Figure 8: Phase diagram of the holographic QCD as seen by [T7]. The figure demonstrates flavor brane 
configurations in different phases. 

density of baryons is exponentially suppressed due to their large mass (proportional to XN C ). This is the 
hadronic phase with no baryons, but with a finite density of mesons and glueballs (lower left corner of the 
diagram in figure [5]). 

At high temperature holographic QCD is in a deconfined phase with restored chiral symmetry. Chiral 
symmetry restoration is realized through a separation of D8- and anti-D8-branes, which is favorable for large 
enough temperature and both zero |30j and finite chemical potential |15j . For moderate chemical potential 
critical temperature of chiral symmetry restoring transition is in general a function of temperature, while 
the critical temperature of the deconfining transition is temperature independent in both conventional large 
N c and holographic models. Therefore there might exist an intermediate phase, in which the matter is 
deconfined, but the symmetry is broken. This indeed happens for sufficiently small values of asymptotic 
separation L between the D8-branes as in the case shown in figure |5] Both deconfining and chiral symmetry 
restoring phase transitions are first order. As noticed by the authors of [17] . in the situation presented in 
figure [5] the shape of the chiral symmetry restoring curve has a peculiar form: it first dips for fi < ms; 
it reaches a minimum for some fi > ms and then raises again. In other words the chiral condensate first 
decreases when \x is increased and then increases again, which is perhaps a similar behavior to some limits 
of QCD, e.g. figure Ha). 

For the values of chemical potential fi > raj in the chiral symmetry broken phase there is a non-zero 
baryon density. As we will explain in the next section this phase should be some kind of quark matter, 
since the density is always large enough for baryons to overlap. If the baryon density is approximated by 
a uniform distribution as in [17] . there is a direct transition from no-baryon hadronic matter to the quark 
matter. The transition would be second order with zero critical baryon density, as shown by figure [3th). 
This is consistent with the observation that the net baryon-baryon force is repulsive in this model. At zero 
temperature the critical chemical potential /i c = ms , while if the temperature is raised the value of /i c will 
slightly decrease. 
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In the finite baryon density phase the geometry of the D8 embedding is modified as described by |17j : 
the brane profile changes from the U-shape to a V-shape as illustrated by figure [8] The new geometry has 
a cusp at the lowest point of the profile. The reason for this modification is the coupling of the baryons 
to scalar fields, mentioned in the previous section. Fluctuations of the scalar fields parameterize transverse 
displacements of the D8-brane profile. In the antipodal case however, only U-shape configuration is possible. 
Indeed in this case the coupling of the scalars to the baryons is forbidden by an additional symmetry and 
thus vanishes, no matter what the baryon density is. 

Studying finite density baryon systems along the lines sketched in the previous section is challenging. For 
generic flat space multi-instanton solution one has to resort to the ADHM construction and solve infinite 
dimensional (in principle) matrix equations, not to mention instanton solutions in curved space, which are 
not known so far. Bergman et al. approximated instanton density by a uniform distribution in three spatial 
dimensions and a delta-function localization at the bottom of the profile z = 0. However as we see in the 
analysis of the previous section, baryon positions are determined by the gravity potential, which prefers to 
have them at z = 0, and the Coulomb-like repulsion of individual baryons. If the baryons are very densely 
packed, it is conceivable that the Coulomb repulsion will overcome the gravity force and push the baryon 
density out to finite values of z. In this respect a more appropriate approximation was considered by Rozali 
et al. In the work |18j they also considered instanton density as uniform in spatial dimensions, but the 
distribution in the holographic dimension was determined dynamically. A non-trivial distribution of density 
was discovered, with a sharp edge at some finite value of z. It was shown that the expansion into the 
holographic direction goes along with the increase of chemical potential. Rozali et al have interpreted the 
finite density in the holographic direction as the Fermi sea with the sharp edge - Fermi surface. The Fermi 
sea must be a Fermi sea of quarks, because in this description baryons are neither fermions, nor bosons. 

Although uniform density is an approximation it is likely an appropriate description of the quark phases, 
when the densities are high. One can expect that in such a case we deal with a weakly interacting quark 
fluid. For smaller densities, such that the size of the baryon a is much smaller than the average inter-baryon 
distance, a better approximation may be required. Recall that in the large N c limit and chemical potential 
fi ~ tub baryons should form a crystal, see explanation in section 12.21 Therefore a better description 
at small densities should be provided by periodic arrays (lattices) of instantons. We will provide such a 
description in terms of the generalized Sakai-Sugimoto model summarized in section [3] and investigate the 
phase space in the intermediate regime of small densities. We will probe a small window of values of the 
chemical potential around /i = ms, where very interesting physics resides. In particular, in this window 
we can observe transitions to structures with finite width in the holographic dimension, which we interpret 
as an onset of the quarkyonic phase (the terminology was explained in section |2~21 . In this respect we will 
supplement the phase diagram in figure [8] by zooming on the values of /i ~ tub ■ 

We also expect that by looking at baryon densities for which it make sense to talk about lattices we 
should be able to establish connections to the results on finite density Skyrme model reviewed in section [2. 31 
Recall that skyrmion lattices exhibit transitions to the new ones with half-skyrmion symmetries. Since there 
is a natural relation between holographic baryons and skyrmions, one can expect that the features of the 
skyrmion lattices should be inherited by the lattices of holographic baryons. So what are the skyrmion-half 
skyrmion transitions in terms of instantons? To answer this question in principle one can use the same 
approach as used below, but this will involve solving for a 3D instanton lattices, which is beyond the scope 
of the present work. A qualitative explanation was proposed by Rho, Sin and Zahed in |20j . Let us briefly 
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summarize the idea. 

As will be clear from section [6] a periodic ID chain of instantons on R 4 with the same orientation twist 
between all nearest neighbors, or equivalently a single instanton on R 3 x S 1 with a non-trivial holonomy 
around S 1 is similar to a pair of 3D BPS monopoles if the instanton separation (radius of S 1 ) is small 
enough [55]. For the twist angle <fi = ir the monopoles have equal mass and opposite electric and magnetic 
charges. The separation between monopoles depends on the radius of S 1 . When the radius is small (high 
density) the topological density is independent from the coordinate on S , but localized in two well separated 
points (monopole centers) in R 3 . For the values of the S 1 radius comparable to the instanton size and larger 
the topological density is localized at one point and does not look as 3D solution any more. 

The authors of [20] have conjectured that the high density splitting of instantons into monopoles is akin 
to the skyrmion-half-skyrmion transition. This conjecture is qualitatively supported by numerical studies of 
skyrmion chains in |59j : at large densities the chain becomes thicker in the transverse directions exhibiting 
skyrmion constituents with fractional topological charge. Some details however remain unclear. In particular 
for a ID chain the instanton-monopole transition is a smooth crossover, while for skyrmions it is expected to 
be a second order transition. We will also demonstrate below that the monopole separation does not depend 
on the baryon density in the same way as for half-skyrmions. Perhaps these issues will be naturally resolved 
in the case of 3D instanton lattices, but the latter have yet to be analyzed. 

4.2 Lattices of holographic baryons 

Let us now explain in more detail why we are interested in the small range of values of chemical potential 
around the critical value /i = tub and what we expect to find there. As we briefly mentioned in the end 
of section 12.21 it is important that we are working in the regime of large A. In this regime the mass of the 
baryon is 0(A) while the baryon interactions are 0(1). One consequence of this fact is that the natural scale 
for the baryon physics is proportional to 1/y/X in units of 1/Mkk- For example, such is the scaling of the 
baryon radius as we saw in section I3~2l As we are going to see later 1/y/X is also a characteristic scale for the 
equilibrium lattice spacing in the holographic crystal. At this scale the baryon lattice undergo a sequence of 
phase transitions. 

In the meantime baryon density pertinent to the scale 1/vA correspond to values of the chemical potential 
in the small vicinity of /i c : 

li = li c + 0(A°) . (48) 

To probe the values of the chemical potential at finite distance away from /i c one has to consider lattices 
with spacing much smaller than 1/vA, that is smaller than the size of the baryon. Accordingly the baryons 
will overlap so that the constituent quarks will not know which baryon they belong to. So this phase must 
be some kind of a weakly interacting quark fluid. In such a case the exact instanton description of the 
baryons becomes tricky, while uniform density may be a decent approximation. In summary nuclear matter 
phase is contained in an asymptotically small window of the chemical potential Afj, ~ jtt c /A, cf. figure [5] As 
seen by the phase diagram on figure [8J this phase immediately changes to a quark matter phase once the 
chemical potential is increased beyond /i c . Nevertheless nuclear matter phase is quite interesting and have 
a rich structure as demonstrated by figure [9] 

At arbitrary low densities we have a normal 3D crystal sitting at the bottom of the potential well in 
the holographic dimension. When the density is increased the baryons will start feeling tight in 3D due 
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Figure 9: When zooming on the small vicinity of the critical chemical potential the phase diagram on 
figure [8] exhibits a richer structure. After the chemical potential reaches critical value fi c baryons are created 
and form a crystal. The crystal becomes multi-layered in the holographic dimension, when the chemical 
potential is increased, where layer production must correspond to the emergence of the quark Fermi sea of 
the quarkyonic phase conjectured in [5]. 



to Coulomb repulsion. It is thus conceivable that at some large density the baryons will pop out from 
the 3D alignment into the holographic dimension, forming a 4D structure. Analogous effect has already 
been observed by Rozali et al. [18j in the approximation of the uniform baryon density. In the latter study 
the 4D thickness of the baryon substance was increasing with the chemical potential and was interpreted 
as a development of the quark Fermi sea. The uniform density of baryons however is not an appropriate 
description for the nuclear matter phase. Instead we need to consider honest lattices of baryons described 
by periodic instanton solutions. 

Working with instantons we will follow the procedure outlined in section 13.21 We will consider periodic 
arrays of instantons for which equilibrium density is related to the spacing d between instantons. We will 
describe the baryon lattice by constructing an appropriate instanton solution in flat space, which corresponds 
to the zeroth order approximation of non-interacting baryons. The zeroth order solution will then be used to 
compute corrections to the free energy due to interactions between instantons. Flat space solutions depend 
on various moduli including instanton positions, size and orientation. Instanton interactions will lift the 
degeneracy. By varying the free energy with respect to moduli we will find their equilibrium values. This 
way we will see that at large enough density the 4D lattice will replace the 3D one in a series of phase 
transitions. 

Although our ultimate goal is to study realistic 3D instanton lattices, in this paper only a ID chain will 
receive the full treatment. For the 3D lattices we will restrict ourselves to the point charge approximation. 
As we will see this is enough to demonstrate the main feature of the holographic nuclear matter: at some 
critical density the 3D baryons will pop out to the holographic dimension. More precisely the lattice will 
split into layers separated along this dimension. Finite instanton size will have the following consequences. 
First of all finite size will reduce the Coulomb interaction of closely packed instantons: the interaction will 
be partially screened. Second, phase space of lattices of finite size instantons have a much reacher structure 
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due to orientation and size dependence of interactions, while the zero-size instantons do not care about their 
orientation. 

In general we expect the following to happen with the 3D or ID lattice of instantons when we vary the 
density (figure [9]) • At small densities the gravity force (curvature) in the transverse dimensions manages to 
keep instantons aligned in the space dimensions despite their electrostatic repulsion. At a critical density 
instanton repulsion will push the instantons out to a direction transverse to the lattice (holographic dimen- 
sion). The ID lattice will split forming a zigzag- like configuration, when every pair of nearest neighbors 
shift in the opposite direction, while the 3D lattice will analogously split into two layers. The transition will 
be second order, with the order parameter being the amplitude of the displacement of nearest neighbors. 
The displacement will grow with density further separating layers in the transverse dimension. For non- 
infinitesimal separation between the layers one can expect other kind of phase transitions between different 
types of lattices related to the change of the dimensionality of the lattices from 3D (ID) to 4D (2D). 

For larger densities (layer separation) the nearest neighbors in each sublattice will become close. One 
can then expect a new transition to a multi-layer lattice (not necessarily to a 3- layer). The transition will 
most likely be first order. Therefore while the density is increased the lattice undergoes a sequence of phase 
transitions changing the lattice structure and multiplying layers in the transverse dimension. When the 
lattice becomes mesoscopic in this direction the instantons will start overlapping in the 3D projection. Here 
the story should match the unform density analysis of |17[ 118]. The baryons will melt into a quark liquid. 
The finite width of the liquid in the holographic dimension will be analogous to the Fermi sea of quarks. 
Therefore we conclude that beyond the first transition in the sequence the phase is resembling the quarkyonic 
one conjectured by McLerran and Pisarski in [5]: a Fermi sea with both quark and baryonic manifestation. 
The first transition just marks the onset of the quarkyonic phase. 

In the remainder of this paper we will consider two toy-models. The first one is a lattice of point charges 
to demonstrate the first transition from a single-layer 3D lattice to a multi-layer lattice in 4D. The second 
toy-model is a ID chain of instantons. To keep such a chain stable we introduce an ad hoc curvature in the 
transverse spatial directions (j44|) so that the chain can only zigzag into the holographic direction. This also 
allows to keep the size of the instanton small a <C d in a controllable way. The ID chain will be shown to 
exhibit the zigzag transition as well as a first order transition to another lattice structure. 

5 First toy model. 3D lattice of point charges 

In this section we will consider an approximation of instantons by point charges. This will allow us to 
demonstrate the instability of 3D lattices against splitting in the holographic dimension. Point charges give 
a good qualitative explanation of what happens with baryon crystals. Taking into account full instanton 
solutions leads to a partial screening of effective interaction between individual instantons and to a richer 
phase space as we will see in the following section, where a ID chain of finite size instantons will be analyzed. 
For a warm-up and for the sake of comparison with the full ID instanton solution in this section we first 
consider the case of a ID lattice of point charges and afterwards proceed to a discussion of 3D lattices. 
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Figure 10: Splitting of the ID chain into the zigzag configuration. At low densities the chain of baryons 
is straight. When the density is increased the chain splits, forming a zigzag configuration. In the case of 
anisotropic curvature (|45j) the deformation of the chain occurs along the holographic direction. 



5.1 Warm-up exercise 

Consider a ID lattice (periodic chain) of point charges. We can think of point charges as of a zero-size limit 
of instantons. To cook-up a ID lattice in the holographic setup one can turn on potential in all transverse 
directions. Here we will consider a case of anisotropic curvature, i.e. use the coupling (|45|) in the expression 
for the energy (|38[) . In such a case dislocations in the holographic direction z are preferable as long as 
M' ^> M . Recall that we label the direction along the chain as X4 and the transverse spatial coordinates as 
x\ and X2, while X3 = z is the holographic coordinate . The density in this case is 1/d, where d is the lattice 
spacing. 

Due to the gravity force the chain must be straight at arbitrary low densities. When the density is 
increased the instantons will be driven out of the ID alignment as long as resisting to the Coulomb repulsion 
will cost more energy than climbing up the gravitational well. Apparently each pair of nearest neighbors will 
prefer to displace in the opposite directions away from the chain. As a result the ID lattice will split into 
two sublattices separated in the holographic dimension forming a "zigzag" as in figure [TU] From the point of 
view of the D8-branes neighboring instantons move to the opposite branches of the U-shape configuration. 

Let us now study this transition quantitatively. We replace the instanton density I{x) by the sum of 
delta-functions 

00 

J(x) = £ *(0( x — nd) , (49) 



n— — 00 



where d is a 4- vector generating translations from one site of the chain to a neighboring one, here chosen to 
be along X4. For the straight chain the non-abelian part of the energy, calculated per instanton, gives 



E NA = N C XM J dx 3 J d 3 x I{x) (1 + M' 2 {x\ + x\) + M 2 xf) = 


= N C XM (l + M' 2 {x\ + x\f + M 2 xf \ . (50) 
The energy is minimized for the choice x\ — X2 — £3 = 0. The solution to equation (|40[) in the case of the 
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point charges is the 4D Coulomb potential. The abelian energy per instanton is then given by the sum 



C 4AM V (nd) 2 AM I2d 2 ' 1 ' 

Let us investigate stability of the charges in the chain against the transition to the zigzag. The condition 
x\ = X2 = holds since we have frozen these directions by requiring M' M . In the zigzag phase every 
even charge is displaced by amount e in one direction along x% , while every odd one is displaced by the same 
amount in the opposite direction. As a result the average non-abelian energy (per instanton) will be shifted 
x\ -> e 2 : 

£na -> E NA + N c XM 3 e 2 . (52) 
It is also straightforward to evaluate the Coulomb energy per instanton in this case 

Nc ( 1 1 \ _ N c { ir 2 ir 7re\ 

Sc -4AM 2. P)2 + If (nd) 2 + (2e) 2 ~ AM [iSd? + T&d ' (53) 

^cvon n#0 V ' odd n K ' y ' J v 7 

The expression for the total energy can be expanded at small e: 

+ +0W ). (54) 

The sign of the e 2 -term depends on the density (lattice spacing). It is positive for small densities and negative 
for large. Therefore there is a second order phase transition, in which e acquires a non-trivial vev, i.e. the 
straight chain transforms to a zigzag. The critical density corresponds to the point, where e 2 -term changes 
sign: 

d = d c = ^ = . (55) 

2-3 1 /4AiVA 

For the spacing slightly smaller than d c the order (zigzag) parameter has the mean field behavior: 

(e) ~ ±^ y/d c (d c - d) . (56) 
More generally e is given by the solution to 

(7re/d) 3 cosh 2 (7re/d) = 3/4\ 4 , , 

sinh(27re/d) - (27re/d) 4 \ d ) ' ( ' 

Graphical solution to this equation is presented in figure[TTJ The zigzag amplitude grows with density. Notice 
that the units for the critical spacing d c are 1/(M V / A) - the same as for the equilibrium size of instantons 
derived in section [3~2| 

One can further study the fate of the zigzag in this toy-model. Once the zigzag amplitude is large enough 
one can think of the zigzag as of a 2D lattice with 2 layers in the holographic dimension. For larger densities 
the number of layers will grow, although the sequence of transitions is not necessarily trivial. More precisely 
if we study a sequence of transitions up to 3-layer lattices we will find that the straight chain turns first to 
the zigzag and then to 3 layers. However if 4-layer lattice configurations are also included in consideration 
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d c 

Figure 11: Zigzag amplitude e as a function of the lattice spacing d for point-like instantons. 
the sequence becomes as follows 

1^2^4^3^4, (58) 

where numbers indicate the number of layers in the lattice and the corresponding configurations are sum- 
marized in table [TJ 

Table [1] demonstrates configurations that appear in the sequence of transitions of the ID chain in the 
case one restricts to the lattices with no more than 4 layers. It shows a formula for the free energy (per 
charge) and a critical density at which the transition to this configuration occurs. The critical density is 
shown in units of the density of the zigzag transition, which from (|55|) can be found to be p c ~ 0.84 in units 
of My/\. The apparent 3- layer configuration is shown in the fourth row of table [TJ If we just have restricted 
our attention to configurations with at most 3 layers, then after the zigzag transition we would have a first 
order transition to the 3-layer lattice at p ~ 2.3p c . In the equation for the free energy of 3 layers we denote 
e the distance between nearest layers as shown in the corresponding diagram. It turns out however that 
the transition to 4 layers occurs at smaller densities, p — 2.26p c . Interestingly the free energy of the 4-layer 
lattice is smaller than that of the 3-layer one for 2.26p c > p > 2Ap c and for p > 3.08p c , while in between 3 
layers win. In both cases the 4-layer configuration is the same, as one can see from the third and the fourth 
rows of table [TJ although in the former case it looks more like a deformation of 2-layer lattice, rather than 
a 4-layer one. That is why in the formula for the free energy we choose e to be a distance from the central 
axis to the imaginary 2 layers, while the parameter 5 measure the amount of deformation, or the "zigzag 
parameter", of imaginary layers. All the transitions in table[TJare first order, except for the zigzag transition, 
which is second order. 

In figure [T^J we plot the free energies of various configurations as functions of the density and summarize 
the phase diagram of the ID chain of point-like instantons (with the restriction of 4 layers). As can be inferred 
from the left part of the figure the energies of various configurations are often very close. For example, for 
p c < p < 2.4p c , the energies of the zigzag, 3-layer (and also 2-layer) configurations are impressively close to 
each other. Although we have not studied the transitions beyond 4 layers it is clear that the sequence of 
transitions will not be trivial. Specifically it looks likely that the transition to 5 layers will occur before the 
second transition to 4-layers. Complexity will grow fast with the number of layers and present a challenge 



34 



Configuration 



Free energy 



x 3 



iV c ^ 2 p 2 



AM 12 



x 3 



^AMV+J^ f^ + ^tanhTrep 



16AM V 3 e 



2.26p c 




N C \M 3 (e 2 + (S 2 ) + 



iV fi 



32AM V 6 



7T 2 p 2 



irp irep 7rptanh7r(e + S)p/2 

coth — — + — h 

e 2 2(e + S) 

7rp , nSp 7rptanh7r(e — 8)p/2\ 

— tanh — - H ; v „ x ;p/ 

<5 2 2(e-<$) / 



x 3 * 



2.4p c 



■ N r XM 



3^2 



36AM 



tt 2 /? 2 



A'4 



7rp , 2nep Attp , 7rep 
— coth — — - H tanh —f- 



x 3 L 

e 



3.08p c 



same as in the 3 row 



-e+60- • - • • • • • 



X4 



-e-<5 



Table 1: Various n-layer conhgurations that appear in the sequence of transitions in the ID lattice of point- 
like instantons. The 2 column and the 4 — column contain the critical density p c and the free energy (per 
instanton) respectively for a configuration shown in the 3~ column. Critical density is measured in units of 
the density of the zigzag transition given by (|55"|) . Notice that lattices in the 3~ and 5~ rows are the same, 
albeit at different densities. Only conhgurations with no more than 4 layers are analyzed. 
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Figure 12: (a) Free energy (per instanton) of various ID point charge configurations as functions of the 
density in units of M\X. (b) Zero temperature phase diagram of point-like-instanton chain. Dashed line 
indicates the second order phase transition, while the solid lines are first order transitions. 

for numerical analysis. 

This otherwise trivial ID exercise demonstrates main features of the crystal of holographic baryons. When 
squeezed enough the baryons cannot seat in the regular space dimensions and pop into the holographic one. 
The phase space of the crystal has a very rich structure as it involves various crystal configurations in the 
additional dimension. In section |6] we promote the chain of point charges to the full instanton solution. 
Interference between finite size instantons will effectively reduce their repulsion and therefore increase the 
critical density. The phase space structure will also depend on the mutual orientation of instantons. But 
first, let us slightly elaborate on the point charge case and consider a 3D lattice. 



5.2 3D lattices 

Let us start our discussion from the simple cubic (sc) lattice. Similarly to the ID lattice we expect that at 
large density Coulomb repulsion will expel instantons from the 3D alignment out to the holographic dimen- 
sion. The lattice will become 4D with multiple number of 3D layers. The most probable first transition must 
be the one, in which any pair of nearest neighbors displace in the opposite directions along the holographic 
dimension, in analogy with the zigzag transition. For the sc lattice this transition corresponds to every even 
site moving one way and every odd one in the opposite. In this second order phase transition the original sc 
lattice will split into two fee layers. 

The above transition can be studied quantitatively. The analysis of stability proceeds similarly to the ID 
case. The only difference is that in 2 and 3 dimensions the Coulomb energy Ec per instanton will diverge. 
In the case of a 3D lattice it is linearly divergent: 

£ dHnl+nl + ni)* 1 -' L "> °° ' (59) 

(ni,7l3,Tls)^0 v 1 1 i ' 

where L is the linear size of the system. However for the question of stability one does not need to compute 
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the sum itself, but rather its variation in case of the lattice expansion to the 4th dimension. In other words 
we are not interested in the infinite constant Eq, cf. (|54p . but rather in the finite e 2 and e 4 terms. It is 
easy to show that to find the energy shift for a split configuration one needs to compute the following two 
quantities 

/ j / 9 9 9 \ 2 ' ^ / j / 9 9 o \ 3 ' \^^/ 

oddK+ n 2 + «3) odd (ni+ni + ni) 

where the sums arc taken over odd sites: ri\ + n% + H3 - odd. Numerical evaluation of the sums gives values 
A/i 2 ~ 10.0, t ~ 6.60. The regularized expression for the energy per instanton of the split lattice is 

The critical density corresponds to the lattice spacing 

0.69—^, (62) 
while the generalization of equation (l56l) for the absolute value of the order parameter reads 



< c )~±-|Vdc(dc-d). (63) 

Thus for a sc lattice with the lattice spacing smaller than (|52l it is favorable to split. In reality however 
the sc phase of repelling instantons is not the most favorable configuration at a given density. The minimum 
energy configuration must correspond to close packing. Close packing gives the largest interatomic distance 
between the nearest neighbors for a given density. In the 3D case this is achieved in the fee lattice. 

Stability analysis of the fee lattice is trickier because for the fee lattice there is no natural way to split 
into two sublattices. One can expect two things to happen when the density is increased. First there can 
be a first order transition to a multi-layer 4D lattice. Second the splitting into two sublattices can occur 
through a breaking of the cubic symmetries (half of nearest neighbors will go one way and another half 
the other way). The latter will be restored as soon as the separation between the sublattices will be large 
enough. Let us also mention the bee lattice. For the latter there is a natural splitting to two sc sublattices. 

We expect the finite density phase space of 3D lattices to be qualitatively similar to that of ID lattices. 
When the lattices are squeezed the instantons will pile-up in the holographic dimension producing more 
layers and rearranging the lattice structure. In the meantime the former phase space will be richer as more 
lattice configurations are possible. For small thickness of a lattice in holographic dimension its structure 
will be defined by the 3D intuition, but when the number of layers is sufficiently large the lattice will start 
exhibiting its 4D nature. 

In this section we have shown that lattices of point-like instantons (baryons) undergo popcorn transitions 
to configurations with finite thickness in the holographic dimension. In other words the boundary theory 
lattices become bulk lattices. The characteristic scale for such transitions is provided by the same combination 
l/(s/XM) as in the case of the baryon size. In fact the same scale will remain for the transitions in the case 
of finite size instantons as we are going to see in the next section. This is an indication that the transitions 
occur when the baryons start overlapping. We will continue this discussion in the concluding section [7] after 
the analysis of the effects of finite instanton size. 
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6 ID chain of instantons — a second toy model 



In this section we will analyze a model of a ID crystal (chain) of baryons. In contrast with the previous 
section the full instanton solution will be used for the chain. Full solution introduces additional parameters, 
such as instanton size and orientation. We start from a straight chain: this solution was originally found by 
Kraan and van Baal [50] ; and compute corrections to the energy of the flat space instanton. The corrections 
stabilize the size of the instantons and their orientation angle. In particular we show that the preferable 
configuration at low densities is "antiferromagnetic" . In the case of a generic instanton configuration the 
Coulomb energy integral cannot be computed analytically. Instead of evaluating it numerically we restrict 
our analysis to the case of small instantons, where analytical derivation is possible. 

Next we derive the zigzag solution to the ADHM equations. We show that at large densities the zigzag 
configuration is preferable over the straight chain(f| At small values of the zigzag amplitude the twist angle 
remains tf> = tt. By analyzing more general twist configuration (with non-abelian twist) we show that at 
larger zigzag amplitude the twist angle becomes (f> = tt/2. 

6.1 Straight chain 

6.1.1 Straight periodic twisted chain of instantons 

A solution for an infinite periodic chain of SU (2) instantons was first obtained by Harrington and Shepard 
in |61j in the context of finite temperature field theory. This solution corresponds to parallel-oriented 
instantons. It was later noticed by Rossi [62 that through a gauge transformation it can be transformed into 
a single 3D BPS monopole solution. In particular the instanton density does not depend on the coordinate 
along the chain. This instanton-monopole equivalence was generalized by Lee an Yi through a study of 
moduli space of calorons - instantons on R 3 x S 1 [58]. More generally a SU(N) caloron solution with a 
non-trivial holonomy around the S 1 is equivalent to ./V BPS monopoles. Net electric and magnetic charges 
of the monopoles vanish and the net topological charge (mass) is that of the instanton. The holonomy of an 
instanton on i? 3 x S 1 has an interpretation of a constant relative orientation twist between the instantons 
in the infinite chain in i? 4 . An explicit solution for a twisted instanton chain was found by Kraan and van 
Baal in [60]. Let us briefly review this result. 

To find a solution for parallel instantons it is enough to write down the 't Hooft ansatz. In the case of 
instantons of variable orientation one has to resort to the details of the ADHM construction [49] . ADHM 
data are matrices, which contain the information about the locations of the instanton centers, their radii and 
orientations. For an infinite chain the matrices are infinite dimensional. More specifically, Sp(k) A-instanton 
solution is encoded by two quaternionic matrices: an N x N symmetric matrix X and a k x N vector Y, 
which satisfy the constraint 

X^X + is real symmetric. (64) 

In the case of Sp(l) ~ SU(2) the solution can be reformulated in terms of four real symmetric matrices 
r M and real vectors Y^: 

X = Y^t^ = T 4 + ir 3 T j , Y = Y»t» = Y 4 + ir^Y 1 , (65) 

8 This statement is of course true only for the densities, which are not too high. At higher densities the zigzag will yield to 
multi-layer structures as in the case of point charges. Remarkably even the straight chain is more favorable than the zigzag at 
asymptotically high densities as we show in appendix Icl 
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where t } , j = 1, 2, 3 are Pauli matrices. The real matrices satisfy the constraint 

([V**, V] + <E>Y U - Y v ®Y») = - e^ ukX ([r K , T A ] + Y K ® Y x - Y x ® F K ) , (66) 

where (F M ® Y"") mn = Y£ x Y£. Physically, the diagonal matrix elements F^ n are the 4D coordinates of the 
instanton centers, Yfi combine the radii and the SU(2) orientations of the instantons, while the off-diagonal 
matrix elements are determined by the condition (|6"6"|) . (Modulo a common O(N) symmetry of all the 

r&„ and y£.) 

For the purposes of this paper there is no need to find the expression for the gauge potential itself. What 
we only need to know is just the expression for the instanton density. In terms of the matrices T' 1 and it 
can be constructed as follows. Provided that (|66f is satisfied define a real symmetric N x N matrix 

L(x) = (x"l -r")(a^l -T") + Y»®Y». (67) 

The instanton density is then given by |63j 

I(x) = --^□□logdet(L(x)). (68) 

We are interested in constructing a solution, which corresponds to an array of equal size instantons 
arranged in a ID periodic lattice. Periodicity means discrete translational symmetry S : a; M — > x^ + d^ of 
the ADHM solution; in the language of the and Y£, this symmetry acts as 

r M ->■ S'^S = + d^l ((to keep the i"l -T^ invariant)) (69) 
(Y*t») Y,G{Y^)S mn = (yy), (70) 

m 

for some 0(N = oo) matrix S mn and SU(2) matrix G. 

Physically, G rotates the orientation of an instanton relative to its immediate neighbor. The more distant 
neighbors are related by G n ~ m rotations, which generate a U(l) subgroup of the SU(2). Without loss of 
generality, we take G = exp(i0T3/2) for some "twist" angle <j> between and 27r. We also take the direction 
of the instanton chain to be the x while the transverse directions are x , x 2 , and a; 3 = z; in terms of 
equation (IM)) this means = (0, 0,0, d). Finally, we take S m , n — 8 m ,n+i (shifts from the n— instanton 
to the (n + 1)—) and Yq — (0,0,0, a) (where a is the radius of the 0— instanton). Consequently, the 
translational symmetry (|69j) -(|70 j) requires 

iT»Y£ = aexpfmfrg) <^ Y£ = (0,0,osin(r^/2),ocos(n^/2)) (71) 
and r^„ = d6^xnS mn + f "(m - n) , (72) 

where the F M (m — n) do not have separate dependence on m and n but only on m — n. Combining these 
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symmetry conditions with the ADHM constraint (1551) . we get 



® Y£ = a 2 cos [(m - n)0/2] , (73) 
r tnn = dnS mn , (74) 

r m „ = r 2 m = o, (75) 

o a 2 sin f(m — n)6/2] „ „ . . 

r mn = -T X for m + n, but for 772 = 77. 76) 

d 777 — 77 

To calculate the instanton density of the periodic chain we need the determinant of the infinite matrix 
L (I67[) . This determinant is badly divergent, but we may obtain it up to an overall infinite-but-constant 
factor from the derivatives 

d M logdet(L(a:)) = 2tr((^l-r M ) J L- 1 (a;)). (77) 

For 3 of these derivatives (/x = 1,2,3) the trace converges while for /i — 4 the trace diverges but can be 
regularized using symmetry — > — 24, 77 — > —n. To evaluate those traces, it is natural to use Fourier trans- 
form from infinite matrices to linear operators acting on periodic functions of 9 (mod 27rj^|. Consequently, 
L becomes 

..2 



L = x{ + xi + (x-3-r^)) + \x± + id— j + T(6), (78) 

where T(fi) = ira 2 5(9-^) + 7ra 2 <5(0+^), (79) 

.1-4- for < 6» < I and 2tt - I < 6» < 2tt, 

and r 3 (6-) = — x ^ 2 " 2 2 (80) 



d 



-f- for § <0<2tt-& 

2tt 2 2 1 



and L 1 becomes the Green's function of this operator. Calculating this Green's function one obtains the 
following expression upon integration of the traces (|77[) : 

dct(L) = f cosh + ^ sinh ^ f cosh ( 2 * ~ ^ + ^ sinh (2 " " ^ 



d dri d J \ d dri d 

r\ + r 2 - (™ 2 /d) 2 ^ 0n ^ (2tt - 0)r 2 



2nr2 

27TX4 



sinh — — sinh ■ 



where 



2d 

This is precisely the result obtain by Kraan and van Baal in 

9 This operation is also known as Nahm transform |64l . 



(81) 



a^-270 
3+ 2d 

*3 + ^y. (82) 
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6.1.2 Total energy of the straight baryon chain 

Let us now apply the above instanton solution to the holographic model of baryons. The goal is to find the 
energy of the corresponding multi-baryon configuration. As outlined in section [3] to find the energy up to 
O(X ) order we need to plug the O(X) (flat space) instanton solution of interest into equation ([55)1 and then 
minimize the energy with respect to the moduli. To stabilize the chain in ID we introduce curvature in all 
directions transverse to the chain. 

The expression for the instanton density follows from (|68|) after plugging solution (|81j) . The non-abelian 
part of the energy (|38j) reduces to a calculation of moments of the instanton density. Although the instanton 
density following from (|81l) is way too complicated to print, the integrals 

(x 2 ) = J^Ax A J A^xxj x I{x) (83) 
can be calculated rather easily by integrating by parts. This gives us 

(4) = <* 2 > 2 = y (84) 

— exactly as for a single stand-alone instanton of radius a — but 

(4) = y + 1^x0(2^-0), (85) 

where the second term is due to interference between the instantons. Curiously, the interference term vanishes 
for <j> = 0, i .e. for instantons with the same SU(2) orientations. 

For a gauge coupling rising isotropically in all directions transverse to the instanton chain, as in the case 
described by equation (jUJ) , and instantons of small size a -C M _1 , equations (|84|) and (|85|) give us the 
non-abelian energy (per instanton) of the Kraan-van Baal chain as 

£ NA = N C XM (l + ~M 2 a 2 + x 0(2tt - 0) + 0{M 4 a 4 )j . (86) 

To compute the Coulomb energy one first has to find a solution to (|4"0)) . In light of equation (|6"8"|) 

A (x) = — — □logdet(£(a;)) + const. (87) 
Consequently, the abelian Coulomb energy per instanton is 

N r XM f d 



E C = ^fJ o dx 4 Jd 3 x(d,A f 

J dx 4 J d 3 .T (df.D logdet(i)) i 



N c <- d 



256tt 2 AM 

For generic lattice spacings d, this integral is too hard to take analytically. But it becomes much simpler 
in the d 3> a limit of well-separated instantons, and also in the opposite d <C a limit, where the instantons 
merge into a continuous line. Here we present the results for the energy in these two limits and refer to 
appendix [A] for more details. 
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For large lattice spacing d^> a the Coulomb energy (|88p evaluates to 

N, 



E ° XM 



2 i <i(„ A\2 

(89) 



1 Air 2 + 3(tt - y; 2 4 

5^ + 30^ + ° {a /d } 



Combining it with the non-abelian energy (|86[) and minimizing with respect to the instanton radius a and 
twist angle <fr, we find the minimum at 

2 3 

0[@min] = 7T, a[@min] = a - ^5 + °( a o/ d2 ) i ( 90 ) 

where ao is the equilibrium radius of a standalone instanton (|46[) . 

In the opposite limit d <C a of densely packed overlapping instantons, the instanton density becomes 
independent of the £4 coordinate, while in the 3 transverse dimensions it becomes concentrated at two 
widely separated points 

= l ' '^ 27 ""^) and * [2) = (9i) 

note that the distance b = ira 2 jd between these points is much larger than the naive instanton radius 
a. Effectively the solution looks like a pair of monopoles on R 3 , in accordance with the analysis of [55]. 
Monopole centers sit at X^- 1 ' and X^ 2 \ The net topological charge of the monopoles is equal to one - 
instanton charge per unit cell, while the distribution of the charge between them depends on the value of <p. 
Indeed for d -C a the 3D monopole density takes the form (see details in appendix |A} 

h D {X) = dxI 4D (X) a ix^'(I-iW) + ^x^(i-i< 2 ») + 0(e~ 2 ^/ d ). (92) 

2n 2tt 

This density holds for well separated monopoles such that their topological charges are <p/2-K and (2ir — 4>)/2tt 
respectively. One can also show that the monopoles have opposite electric and magnetic charges. In [55] 
the monopoles are obtained from an instanton solution on R 3 x S 1 — caloron. By comparing (|92p to the 
monopole charge densities in |58j one concludes that the twist angle <f> must be identified with the holonomy 
of the caloron solution around S 1 . 

The Coulomb energy of the large density solution can be expressed as 

Ec -8^MdW 1 + ^2^ + b )' (93) 

where pi and P2 are effective radii of the monopole density concentrations at X^ and X^ 2 \ Formal 
expressions for them are given by equations (| 1 75[) and (I178P in the appendix [A] Plugging these charge radii 
into equation (|93[) , we write the Coulomb energy of large instantons as 



Here C\ « 1.174 and C2 ~ 1.761 are constants defined by some dimensionless integrals (|176l) and (| 1TT|) . 
Combining this Coulomb energy with the non-abelian energy (|86p and minimizing with respect to the 
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instanton radius a and twist angle </>, we find the minimum at 



7T, 



1/3 



0.842 aZ d 



2/3 ,1/3 



(95) 



Note that when the lattice spacing d becomes much smaller than the equilibrium radius ao of a standalone 
instanton, the formal radius parameter a of the instantons in the chain decreases as a a a^d 1 / 3 . However, 
the actual concentration of the instanton density shrinks to two separate 3D zones of a much smaller size 
O(d) <C a but separated by a much larger distance b — ira 2 /d that increases for small d as 



2.23 a, 



4/3 



dV3 



(96) 



In terms of the root-mean-square (rms) density radius 



<r>, 



d x(x ± + x 2 + x 3 )I 3D (x) 



1/2 



= a x 



4d 2 



1/2 



(97) 



we have something like the blue line in figure [T"3l Here Is,d{x) = d x I^dIx) and exact expressions for the 
moments are used, cf. equations (l84l and (|85l) . Equation (|97|) — and hence the red line on the graph 
- presume <p = 77 f° r au lattice spacings d. For large and small d we have seen that this is indeed the 
equilibrium (i.e. the lowest-energy) value of the twist angle. For the intermediate d ~ ao an argument is 
given in the appendix IBI that </)[@min] — ir, independent of lattice spacing. 

Let us briefly comment on a possible relation between these results and skyrmion-half-skyrmion phase 
transition discussed in section l2~3l Rho et al. [20^ have suggested that the splitting of each instanton into a 
pair of monopoles, which occurs in the large instanton regime d < a, is equivalent to the latter transition. 
Indeed, we have proven that the preferable twist angle is <f> = ir and therefore each monopole carries a half of 



43 



the instanton charge On the other hand the splitting is merely an alternative description of the instanton 
solution related to the fact that for a > d the instanton density is independent from It is a smooth 
crossover rather than a phase transition, and we do not expect chiral symmetry restoration in the sense of 
section 12.31 This is perhaps a drawback of the ID model, which will not be present in 3D lattices. If we 
look at the parameter b of the monopole separation, it smoothly interpolates between a a^/d behavior for 
a <C d and a^ 3 /d 1 / 3 for a 3> d, that is it smoothly grows with density. For the 3D lattice case the smooth 
behavior would break the discrete symmetries. It is then conceivable that in the latter case the crossover 
will be replaced by a phase transition. 

6.2 Abelian zigzag. Exact solution 

Our main interest in the instanton chain is not in its equilibrium configuration but rather its stability (or 
instability) against transverse motion (in £1,2,3 directions, where X3 = z) of the instantons that would get 
them out of the linear alignment. If such motion can decrease the chain's net energy, than the chain is 
unstable. 

In the analysis of stability it will be useful to break the symmetry in the directions transverse to the chain 
and introduce anisotropic gauge coupling (|45l) . The main motivation for introducing the anisotropy is to 
make the chain more unstable in the holographic X3 rather than transverse spatial x\, £2 directions. For this 
apparently we must impose M < M' . The anisotropy also breaks the degeneracy between different directions 
n of the SU{2) twist exp(i0(n • r)/2) between adjacent instantons. The lowest-energy direction of the twist 
is now T3 - which is precisely what we have used in our formulae in the previous section. Specifically, while 
the Coulomb energy does not depend on the twist direction, the non-abelian energy is minimized when the 
largest (a; 2 ) of the instantons is oriented in the lowest-potential direction 2:3. According to equations (|84|) 
(|85p. the instantons are larger in the direction of the twist than in other directions, hence the non-abelian 
energy 

E NA = N c \m(\ + a 2 M' 2 + ^a 2 M 2 + ^^zA x (M l2 {n\ + n\) + M*n%) + 0(a 4 M 4 )) (98) 

is minimized for n = (0, 0, ±1), i. e. SU (2) twist in the r 3 direction. 

Also, M' > M reduces the equilibrium size of standalone or far-apart instantons from (|46p to 

Consequently the anisotropy, apart from stabilizing the chain in the x\, xi directions, provides an additional 
control over the equilibrium size of instantons. This will allow for analytical analysis of phase transitions 
below. 

Before we start the stability analysis we must describe the possible displacement in proper ADHM terms 
of r^ n matrices and Y£ vectors. Moving the instantons' centers in £1,2,3 directions without changing their 

10 The reader should not be confused by the fact that the splitting in our model occurs along the holographic direction. This 
can be cured by a different choice of the twist direction. All twist directions are degenerate in the case of isotropic gauge 
coupling [144 1 . In a 3D lattice the twist direction and thus the instanton splitting, if it occurs, will necessarily be along the 
spatial directions. 
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Xi locations along the chain or any radii or SU (2) orientations means keeping 



Yg = (O,O,asin(n0/2),acos(n< ? !)/2)), T A rnn = dnx6 mn , (100) 
exactly as for the straight chain but changing the T^' 2 ' 3 matrices 
-> r^„ [straight chain] + ST 3 mn , T 1 ^ 

^ ^ ran 5 ^ mn ^ ^-^mn 5 (-^l) 

in a manner that preserves the self-duality equations (|66[) . The solution corresponding to the displacement 
in the X3 direction reads 

ST^l = 0, 5r 3 „„ = 8 mn x <5A 3 [n] . (102) 

The gauge coupling keeps the instanton centers lined up along the X4 axis for low instanton density. At 
high density, such alignment becomes unstable because the abelian Coulomb repulsion between the instantons 
squeeze them out in the directions transverse to the chain. Since the repulsion is strongest between the nearest 
neighbors, the leading instability should have adjacent instantons moved in opposite ways forming a zigzag 
pattern (figure [TU|) . 

SX 3 [n] = e x (-1)". (103) 

After defining the ADHM matrices we make a Fourier transform to map the infinite-dimensional matrices 
to differential operators on a circle. In the following it will be natural to combine wave functions if) (6) and 
ip(9 ± 7r) on the circle into a two-component wave function. Two-component functions provide a natural 
description for a two-layered chain. For some particular cases of the twist angle (f> = or it all expressions 
take particularly compact form. More generally for n layers n-component functions should be used instead. 
Here we choose two-component functions with the following boundary conditions 

(tb(0) + *b(6 + tt)\ tt tt , , x , . . 

• m -(S«-S. + w)J' -i-'-f •<-/»>-*•(-*/». do-.) 

where Ei,E2,E3 are Pauli matrices acting on two components. In the space of such two-component wave 
functions, the T 3 operator becomes 

r 3 = exE 3 + — e(-f<0<f)xE 1 + — ^_e(-f <*<f)Jxi, (105) 

where B is the step- function, i.e. = 1 if — cj)/2 < 6 < cf)/2 and O = otherwise; and we restrict <fi to the 
values < 4> < tt. For cj) — tt the latter expression takes the form 

9 

TTQ, 

r 3 = — x Ei + e x E 3 . (106) 
Matrix T = ® Y^ in the two-component notation becomes 

T(9) = ^ (1 + Ex) (s(6 - |) + 6(0 + |) V < < tt . (107) 



For = 7r the latter can just be written as T = na 8(8 ± 7r/2) x 1, c/. equation ([79 
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Next we plug the data above in the formula for the operator L (|67[) . To calculate the determinant of 
this operator, one needs to calculate its Green's function L^ 1 (6,6q) and consequently the traces f|77j> . As in 
the example of the straight chain the naive determinant is divergent and needs to be regularized, which is 
done by extracting an infinite but constant prefactor. One can reconstruct an analytic expression for det(L) 
with any value of <f>. However only for 4> — and tt the expression is compact enough to present in the 
paper. Accidentally = 7r is the minimal energy configuration of instantons for both the straight chain at 
low density and in a zigzag phase provided that the zigzag amplitude is not too large. A proof of this fact 
will be presented shortly. In the remainder of this section we discuss the case <j) — tt. One finds that the 
regularized determinant reads 

det(L) / irri ira 2 nri \ ( 70*2 7ra 2 71T2 
= cosh — — + — — sinn — — x cosh — — + — — sum 



const \ d dri d J \ d dr 2 d 

r? + r 2 — (7ra 2 /d) 2 , ttti , 7rr 2 2^x4 

+ — — — smh — — smh — — — cos — — 

2r\T2 d d d 

nx 4 ( ( n , tt^i m 2 . , 7rri \ / irr 2 na 2 irr 2 

sin v x cos — — x 2 cosh — — + — — smh — — — 2 cosh — — + — — smh 



where 



d \ \ d dri d J \ d dri d 

o / , 7rri , 7rr9 r? + r 2 . , nri , 7r?'2 \ 

+ sin 2 v x -1 + cosh — - cosh — f- i ^ sinh — - !- sinh — - , (108 

d d 2rir 2 d d J 



r\ 2 = x\ + x\ + {x 3 T \b,) 2 , (109) 



b e = v/4e 2 + (7ra 2 /d) 2 , (110) 

v = arctan — 5-7-r = arcsm— . (Ill) 
na^/d b t 

Naively, we would expect the zigzag deformation to have no effect on the width of the instanton chain 
in the x\ and X2 directions while the width squared in the X3 direction should increase by e 2 . And indeed, 
this is precisely what happens for any instanton radius and lattice spacing: evaluating (I83p by integrating 
by parts, we obtain precisely 

<*?> - <*i> - y' (4) - y + % + ^ (112) 

Consequently, the non-abelian energy of the zigzag is, cf. (|52"|) . 

E NA = N C XM ( (|M 2 + M ) x a 2 + M 2 x + M x e 2 J 

= E NA [e = 0] + N C XM 3 x e 2 . (113) 

As to the Coulomb energy, we will restrict ourselves to the d^$> a regime, where analytical calculation is 
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possible. For small, widely-separated instantons, we approximate 



for r\ ~ a <C d but finite e/d, 



det(i) 



1x2 (, 1 2 2n 71-2 



7ra;4 7rr 2 

cos — ; h cosh — — 

a a 



(rj + 2a 2 r\ - x\) + 0(a H /d 4 ) ) 

2 



na nr2 
- — ; smh — — 

2r 2 d d 



7r 2 a 4 



and consequently obtain 



AM 



1 

5^ 



x smn 

d d 



3tt 2 



80d 2 



80d 2 



tanh(7re/d) 
ire/d 



+ 0(a 4 /d 6 ) 



(114) 



(115) 



Remarkably, the e-dependent part of this Coulomb energy is precisely 5 times smaller than the energy of 
point charges in a similar zigzag formation (|54|) . Mathematically, this fivefold reduction stems from the last 
term in equation (|114l) . which accounts for the interference between the instantons. It is not clear why the 
interference between well-separated small instantons has such a drastic effect on the Coulomb energy of the 
zigzag. Anyhow, the net energy cost of a small zigzag deformation e <C d is 



AE nat = AE 



NA 



AE C = N r \M A x e 2 



240# 



^ 6 6 4 

600d 6 



+ 0(e 4 /d 6 ) 



(116) 



This cost is positive — and hence the straight chain is stable — for all sufficiently large lattice spacings 



d > d r = 



/240 MVA 



(117) 



For smaller lattice spacings d < d c , the energy function ()116|) has a negative coefficient of e 2 but positive 
coefficient of e 4 . Thus, for d < d c the straight chain becomes unstable and there is a second-order phase 
transition to a zigzag configuration. For lattice spacings just below critical, the zigzag parameter e is 



V5 



± x y/d c (d c - d) . 

7T 



(118) 



For smaller lattice spacing (but larger than the instanton size), the zigzag parameter satisfies the same 
equation as in the point charge limit (|57|) with the new d c . Notice again that the new critical spacing is 
precisely \/5 times smaller than the one in the point charge limit. Graphically, the zigzag parameter e as a 
function of the lattice spacing behaves as demonstrated in figure [14l Note that this curve can be trusted only 
for lattice spacings larger than the instanton size since it is based on equation (|115[) for the Coulomb energy 
which presumes rf> o. In units of the equilibrium instanton size (|99p . the critical lattice spacing (|117[) is 



rfcrit 



^ M' 



iM 2 



M 2 



1/4 



1.2a 



(for M' > M). 



(119) 



For a highly anisotropic gauge coupling with M' ^> M, the critical spacing is much larger than the instantons, 
which justifies our approximations. However, for near-isotropic couplings with M' « M, the critical spacing 
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Figure 14: The zigzag amplitude e as a function of the lattice spacing d for small instantons, a <C d. 

is only d cr it ~ 1.32ao, and in this regime we cannot be sure our analysis of the zigzag instability is even 
qualitatively correct. That is, we are not sure that there is a transition from a straight instanton chain to a 
zigzag for M' ps M, never mind any details of such transitions. On the other hand, for the highly-anisotropic 
M' 3> M setups, we can trust that the transition does happen at the critical density (|117|) . and that the 
zigzag amplitude below the transition behaves according to the curve (figure [T4")) . at least for the lattice 
spacings not too much smaller than the critical. 

6.3 General twist. Approximate solution 

The zigzag instanton solution considered in the previous section was derived assuming <f> = tt. When the 
density is increased and the zigzag amplitude grows large enough one expects a series of transitions similar 
to the ones discussed in the case of point charges (table [1] and figure IT2"j) . However the phase space of 
the instanton chain must be richer: e.g. one can expect transitions driven by a change of their mutual 
orientation. Indeed, for the zigzag amplitude e comparable to the chain spacing d, it is natural to expect 
that the orientation of instantons will respect the essentially 2D character of the lattice. In particular twists 
between adjacent instantons may appear along multiple directions in the SU(2) space rather than a single 
U(l) subgroup - a non-abelian zigzag. In this section we consider a more general ansatz for the instanton 
twists and show that while for small amplitudes the <f> = tt twist is indeed preferred, for large ones it 
asymptotes to <fr — n/2, that is an anti-ferromagnetic order inside each of the two zigzag layers. Surprisingly 
non-abelian twist configurations are always disfavored, for any values of density in the d 3> a, M' 3> M 
approximation. 

Expressions for instanton solutions obtained for generic twist are quite monstrous. Even though their 
derivation is straightforward in principle the results can hardly be presented in a compact form. On the other 
hand for the sake of keeping our computations analytical we have to resort to the d^$> a expansion anyway, 
when computing the Coulomb energy. Therefore in this section we will employ this expansion already at the 
level of the instanton solution. 
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6.3.1 General observations 

Let us summarize some observations about the instanton size and energy, based on the examples of single and 
multi-instanton solutions considered above. A single stand-alone instanton of radius a centered at Xig 3 = 
has non-abelian energy (of the SU (2) magnetic fields) 

£na = N C XM x (l + (M' 2 + \ M 2 ) x a 2 ) , (120) 

while the U(l) electric fields (due to Chern-Simons coupling to the instanton density) have 5D Coulomb 
energy 

E ° = m*^- < 121 > 

Minimizing the net energy £?na + Eq with respect to the instanton radius, we find 



1 ,/ 2M' 2 5- 1 / 4 



V XMM' V 5(2Af' 2 + M 2 ) M'r»M \JXMM' 



(122) 



and ^ M • \J^Jm' 2 + Up > M + J\^- . (123) 

XN C A V 5 V 2 M'>m V 5 A 

For an infinite chain (straight or zigzag shaped) of instantons with lattice spacing d much larger than 
the instanton radius a, the non-abelian energy per instanton has general form 

E NA = N C XM (l + (M 12 + i M 2 ) x a 2 + M 2 x e 2 + 

+ M 2 4 x A + M' 2 ^r x B + 0(M' 2 a 6 /d 4 )] , (124) 
d z d l ) 

where e is the zigzag amplitude and A, B are some O(l) functions of the SU(2) twists between instantons 
and also of the e/d ratio, cf. equations (|86p and f| 1 13[) . Likewise, the Coulomb energy per instanton of a 
chain or zigzag has general form 

E ° = m (5. + ? + °<" 2 ^>) ■ ("0 

for some function C of the inter-instanton twists and e/d, cf. (|1 1 5|) . Now let us assume that M' 3> A/ 

and a lattice spacing d ~ (y/\M ) _1 3> a. With these assumptions, we may approximate the net energy as 

E nct = N C XM + N C M' yXMM'a 2 + ^7^2) ( 126 ) 
+ Nr _M ((»»«) x («/<)> + x 8 + ^ x C ) (127, 

+ 0{N c M 2 /M'). (128) 

Minimizing this net energy with respect to the instanton radius we find 

XMM' x a 2 = —= + 0(M/M'), (129) 
v5 
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similar to a standalone instanton, and 

E nct = E B + N C M ((\M 2 d 2 ) x (e/d) 2 + 5 ^*^2) ) + 0{N c M 2 /M') , (130) 

where Eb is the energy of a standalone instanton, c/. equation (|123[) . Note that the lattice spacing d and 
the zigzag amplitude e enter this formula via the ratio e/d and the combination XM 2 d 2 , so the natural scale 
for e and d - at which we might have some phase transition(s) - is 

d,e -i— . (131) 

For M' 3> M this scale is much larger than the equilibrium instanton radius (|122[) . which justifies the 
perturbative expansion of the non-abelian and Coulomb energies in powers of 

a 2 M . , 

# ~ m> « L < 132 > 

This concludes our general analysis, but the devil is in details. To work out the phase structure of the 
instanton zigzag, we need to specify the instanton orientations in SU (2) in terms of some moduli, solve the 
ADHM equations for the non-diagonal matrix elements of the Y l / nn matrices, calculate the instanton number 
density I(X), then use it to obtain the non-abelian and Coulomb energies per instanton. Writing those 
energies in the form (|124l) and (I125[) for some functions A, B, and C of e/d and the moduli of the instanton 
orientations, we plug them into equation (|130|) and minimize the combination 

6.3.2 ADHM data 

As before we are interested in infinite zigzag configurations with periodicity d in X4 direction and amplitude 
e in X3 direction. Therefore instanton positions (diagonal elements X!^ = T£ n ) are given by 

X* = nd, X 3 n = (-l) n e, X\ = Xl = 0, neZ. (134) 

In general we expect e, d — 0(1/ VAM), also we allow for a straight chain with e = 0. We assume that the 
radii and the orientations of the instantons share the geometric symmetries of the zigzag: all instantons have 
the same radius a, while the relative SU(2) orientations between neighboring instantons depend on parity 
of n: 

Y 2n+1 = Ui x Y 2n , Y 2n = U 2 x r 2 „_! , U u U 2 g SU(2). (135) 
In 50(3) terms, U\ and U 2 are rotations through the same angle a but around different axes ni and 112, 

U\ = exp(i(a/2)ni • f) , U 2 = cxp(i(a/2)n2 • f) . (136) 
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We find it convenient to reparameterize these SU (2) matrices in terms of two angles <f> and /? and two unit 
vectors p and q that are perpendicular to each other: 



U 1<2 = exp(T«(/3/2)q • r) x exp(z(0/2)p • f) x exp(^(/V2)q ■ f) (137) 

so that the individual instantons' Y n parameters get a rather symmetric form 

Y n = a x exp(m(0/2)p • f) x exp(i(-l) n (/3/2)q • f). (138) 

The parametrization of (|137l) is redundant, so we limit the range of angles to < < it and < f3 < tt/2. 
But even in this range some redundancies remain: For = the p vector is irrelevant, while for = it we 
have Ui = U2 = ip ■ f regardless of the /3 and q moduli. Also, thanks to the sign redundancy of ADHM data 
Y n — > (±) n Y n , for (3 = tt/2 the exchange 

(0,p,q) o (7r-0,q,p) (139) 

has no physical effect. We should keep these redundancies in mind when we seek the lowest-energy configu- 
rations of the instanton zigzag. 

The instanton centers (|134j) govern the diagonal matrix elements of the quaternion-valued symmetric 
matrix X , which we expand in the basis of r-matrices as in equation (|65[) . The non-diagonal matrix elements 
of X follow from the ADHM equations (|r3|) . We can solve these equations in the matrix form, but we find 
it easier to work in a different basis where the coefficient matrices r M , Y^ <g) Y^ become operators acting 
on two-component wavcfunctions of a periodic variable as we defined in equation (|104[) . In this basis, the 
T mn — Re(Y 7n Y n ) matrix becomes a singular local operator, cf. (|107|) . 



T{0) = — (1 + Si cos (3) x {6(9 + J ) + 6(0 - J ) J . (140) 

Likewise, the imaginary quaternionic matrix Im(YXY n ) = if - t mn becomes a singular local operator if - 1(8) 
where 



t(0) - — 



(p(Ei +cos^) + (p x q)E 3 sin/3) x (-{6(6+^)+i6 

+ q E 2 sin /3 x 



(141) 



As to the matrices, the V is the same differential operator as before 



r 4 = -**| , (142) 

while the 3-vector V becomes (2 x 2)-matrix-valued local function of 8 satisfying ADHM equations in 
3-vector form 

-id^fffi) + T(9)xT(6) + t(9) = 0, (143) 
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subject to boundary conditions 



r(0 = -7r/2) = S 3 T(9 = +tt/2) £ 3 . (144) 

These equations may be solved exactly in terms of an auxiliary variable related to 9 via an elliptic integral, 
or rather two auxiliary variables, one for —0/2 < 9 < cj>/2 and the other for 9 < — 0/2 or 0> 0/2 — the 
solutions are discontinuous at 9 — ±0/2 because of the 5-functions in the t(9) term (| 141[) . However the 
exact solutions are too messy to work with and we will focus on the aCd limit of the T(9). 

Instead of expanding the exact solution in powers of a 2 /d 2 , it is easier to solve the ADHM equations (| 143|) 
perturbatively, starting from zero-order zigzag solution (|102[) considered in the previous section. 

(f ) mn = e5 3 xU-l) n ^ fo(0) = efi 3 xE 3 , (145) 

where n 3 is the unit vector in the X3 direction. The first-order term in T(9) satisfies the linearized equa- 
tions (|143[) , namely 

-id^fx{9) + 2en 3 x [s 3 ,fi(0)] + t{9) = 0. (146) 
The diagonal 11 and 22 matrix elements of the solution to these equations are piecewise-constant functions 

of e, 

(riW ) = lr (p cos/3 + ( p x « S3Sin ^) x (i ®(-t < 9 < 2 } ) ' (147) 

where is the step-function, while the off-diagonal 12 and 21 matrix elements are more complicated: for 
-0/2 < 9 < 0/2 

2 

= ^"(-Psns + v cosh(2(e/d)0) - i(n 3 x v) sinh(2(e/d)0)) , 

where v = —y — y-r ( — p± cosh((7r — 0)e/d) + (ri 3 x qj_) sin/? sinh((7r — cj))e/d)] , (148) 

cosh(7re/a) V / 

while for 0/2 < 9 < vr/2, or -tt/2 < 6» < -0/2 



2 



-i<? 3 sin/3n 3 + v' cosh((e/d)(?r - 2|0|)) \ 
± i(n 3 x v') sinh((e/d)(7r - 2\9\)) J ' 

where v' = —. — ^((n 3 x p) sinh(0e/G?) — q± sin,9cosh(0e/cf) ). (149) 

cosh(7re/a) V / 

In our notations, p 3 and g 3 are the X3 components of the p and q unit vectors while p± and qj^ are their 
components perpendicular to the X3 direction; the ± sign in equation (|149l) is the sign of 9. The next order 
in perturbation theory is rather complicated, so for the present purposes we simply let 

f(9) = en 3 S 3 + fi(0) + 0(a 4 /d 3 ). (150) 

In summary we have derived the ADHM data up to the first non-trivial order in a/d. Now we are going 
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to present the results for the energy of the zigzag configuration as a function of the twist parameters. Details 
of this calculation can be found in the appendix |D] The energy needs to be minimized at given density to 
find the stable zigzag configuration. 



6.3.3 Phase space of the zigzag solution 

To find the minimum energy zigzag configuration of instantons we need to evaluate the energy of the zigzag 
per instanton and read off the coefficients B and C defined by the general form of the net energy (|130|) 
assuming M' 3> M limit. Recall that in this limit the small instanton expansion is justified: a -C d. As 
before, one first computes the determinant of matrix L (|67| . then uses equations (|68[) and (|40|) to find the 
instanton density and gauge field Aq 1 and plugs the result in the formula for the instanton energy (|38l) . 
Expansion of the expression for logdetL for the generalized ADHM data (I140j) and (|147p - (1149|) is presented 
in the appendix [D] Calculation of the net energy per instanton for the optimal instanton radius in the 
M' > M limit gives, cf. equations (fT30|) and ([133]) . 

End = E B + g^y^L x J-(moduli) + 0(N c M 2 /M>), (151) 
where T = 3 - ^SUZ. ^ x (2 - pi cos 2 j3 - (p x q)i sin 2 /?) 

+ x (5 - 4sin 2 /3 + 2pi + 2qi sin 2 f>) 

{ne/a) V / 

+ ' mh((20 ~J )£ ^| x (-4 + 4 sin 2 /? + 2pi - 2qi sin 2 /? 
(ire/ a) cosh{ne/d) V 



+ _° (AAf 2 d 2 ) 2 x (ire/d) 2 . (152) 

7T 4 



To find the phase structure of the instanton zigzag, we now need to minimize this energy with respect to 
the moduli e/d, </>, /3, p, and q. Let's start with the /3, p, and q moduli while keeping the ire/d ratio and 
the (f> angle fixed. Since the energy (|152j) depends linearly on sin 2 /3 and cos 2 /3, its minimum with respect to 
/? must be at either j3 — or j3 = tt/2. For ft = we have an abelian zigzag, the q parameter is irrelevant, 
and the energy dependence on the other moduli has form 

^ = 0) = 3 - M^x[2-pi] + ^P><[5 + 2pi] 

(■jre/d) cosh(TTe/d) L J 7r 4 

This function depends on the p unit vector as T = A + B\t\ where the B coefficient is always positive 
(except for = when p is irrelevant). Consequently, the abelian zigzag always prefer pj_ = — the p 
vector should be parallel to the X 3 axis. For this choice, the abelian zigzag's energy as a function of the 
remaining moduli e/d and 4> is 

80(tt - 4>) tanh(7re/cQ A v sinh((2(/) - 7r)e/d) 

-Aabolian(0,e/d) = 3 + 5x — — 4x — — — 

tt z y-Keld) {ne/a) cosh(7re/a) 

+ _^(AM 2 d 2 ) 2 x (ire/df. (154) 
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Now consider the non-abelian zigzags with (3 = ir/2. For such zigzags, the energy depends on the other 
moduli as 

•F(/3=tt/2) = 3 - 4 ^ 2 "^ x(2-(pxq)l=pl+qi 

' tanh(W(i) x f 1 + 2pi + 2qi) + ; inh /j^-j) £ /j) x - 2q 2 



(7re/d) V ^ x ^7 (7re/d)cosh(7re/(i) 

80 

TT 



+ — (AAf 2 d 2 ) 2 x (ire/d) 2 . (155) 



As a function of the vectors p and q, this energy has form T = A + Bp 2 ± + Cq^, where the coefficients B 
and C can be both positive and negative, so at the minimum the p and q vectors must be either parallel or 
perpendicular to the X 3 axis. Since plq, this gives us three choices: 

(Nl): p||A 3 , qlX 3 , 

jr m = 3 _ ^ - + 3 x tanh(7re/d) _ ^ ^ sinh((20 - n)e/d) 
tt 2 ( 7Te /d) (ire/d) cosh(ire/d) 

80 



+ ^(AAf 2 d 2 ) 2 x (ne/d) 2 , 

TT 

(N2): q||X 3 , pll 3 , 

j- N2 = 3 _ ^ - + 3 x tanh (W^) + 2 x sinh((20-7r)e/d) 

7r 2 i^e/d) ( ne /d) cosh(7re/<i) 

80 

(N3) : pll 3 and q _L A 3 , 

80(tt - 0) tanhfre/rf) sinh((20 - ir)e/d) 

>N3 = o — 5 + 5 X — + X 



+ —t (\M 2 d 2 ) 2 x (ne/df, 



(ire/d) {ne/d) cosh(7re/<i) 

on 

+ —j (XM 2 d 2 ) 2 x (Tre/d) 2 , (156) 

TT 



but the first two choices (Nl) and (N2) are physically equivalent in light of the redundancy (|139[) (which 
exchanges p with q and <fr with tt — <j) for /3 = 7r/2), so we shall focus on (Nl) and (N3) only and drop (N2) 
from further consideration. 

At this point, we have three possibilities — the abelian zigzag, and the non-abelian zigzags (Nl) and 
(N3) — so let's compare their energies at some <j>\ < tt/2 and <p2 — ^ — <f>i > tt/2. It is easy to see that 

-^abelian (02 ) < -^3) (02 ) = F^i^l) < ^abelian (<Pl )) (157) 

so the abelian zigzag with <p = (j>2 > tt/2 always wins over the (N3) zigzag or the abelian zigzag with 
4> = 4>\ < tt/2. The only exception to this rule is <f> = tt/2 for which the (N3) and the abelian zigzag 
have equal energy, but that point is not the global minimum of energy — or even a local minimum since 
d J- abelian / d(f> < at <fi = 7r/2. As to the non-abelian zigzag (Nl) it has lower energy at 4>2 > tt/2 than at 
4>2 = tt — 4>i < 7r/2, but for any particular tfi > tt/2 the choice between the (Nl) and the abelian zigzag 
depends on the e/d ratio. 

Thus far, we have narrowed the choice of zigzags to the abelian and the (Nl). Their energies depend on 
the 4> angle in exactly the same way, up to an overall factor of 2. Indeed, we may write the T functions for 
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the two zigzags as 



•^abelianC^) e/d] d) 

F{m){<t>,z/d\ d) 
where Q(e/d;d) 

and H((/),e/d) 

note J-"abelian ~ -^"(Nl) 



= g(e/d;d) + 2Ti{<P,e/d), 
= g(e/d;d) + H(<p,e/d), 
tanh(7re/d) 80 
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independent of <f> , (160) 



+ 2 x 



tanh(7re/d) 



- 2 x 



sinh((20-7r)e/d) 
(ne/d) cosh(7re/(i) ' 



Ti^e/d). 



(161) 
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For a fixed e/d ratio, both zigzags favor the same angle <p that minimizes the T-L function. Solving the 
minimum equation 

(163) 



7T dU _ 20-tt cosh((20 - w)e/d) 



it cosh(7re/d) 
gives us </) as a function of the zigzag ratio e/d shown in figure [TBT a 





1 e/d 

(a) (b) 

Figure 15: (a) Twist angle as a function of the e/d ratio for an abelian or non-abelian (Nl) zigzags, (b) T-L 
at optimal twist angle 4> as a function of the zigzag ratio e/d. 

For small zigzag ratios e/d, the preferred angle is <fi = it for which H = and both zigzags have the 
same energy. Indeed, at that point the two zigzags become indistinguishable since for <j) = tt the /3 angle 
is irrelevant. However, for larger zigzag ratios the minimum moves to (f> < tt and the two zigzags become 
distinct; for very large e/d, the minimum asymptotes to tt/2. The transition between <f> = tt and <f> < tt 
regimes occurs at (e/d) ~ 0.382 where the second derivative 



d 2 n 



8 / 

— 5-fl — (TTe/d) tanh(7re/<i) 

IT 1 V 



(164) 



changes its sign. 

To find which of the two zigzags is preferable in the regime with an optimal angle <f> < tt we plot 
in figure fT5T b) the difference A J 7 = H (|162j) as a function of e/d. We see that T-L is either zero (when 
the minimum is at <fi = tt), or negative (when the minimum is at <f> < tt), but never becomes positive. 
Consequently, whenever the (Nl) and the abelian zigzags become distinct, the abelian zigzag has lower 
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Figure 16: The Phase diagram of the zigzag configuration. Three curves correspond to the free energy 
per instanton (magenta), zigzag amplitude (blue) and twist angle (yellow) as functions of density. The 
plot exhibits a second order (straight to zigzag) and a first order (orientation flip) phase transitions. The 
second order transition is only distinguishable in the behavior of e, which becomes non-zero above some 
critical density. The plot of the free energy has a "butterfly" typical for a first order transition with stable, 
metastable and unstable parts. The energy is measured with reference to the stand-alone instanton energy 
E B . 

energy. 

We summarize the results of this section in figure IT6l In this figure we plot the free energy (per instanton) 
of the chain configuration, zigzag amplitude e and twist angle <j> as functions of the density in units of y/XM. 
At low densities the chain of instantons is straight (e = 0) and the twist between adjacent instantons is 
4> = 180°. At density dr 1 ~ 1.25 the straight configuration changes to the zigzag in a second order phase 
transition. The twist angle remains 180° in this transition. At density d^ 1 ~ 1.48 the relative twist in the 
zigzag changes from <fr — 180° to <j> ~ 117° in an apparent first order transitions. The zigzag amplitude also 
jumps across the transition. When the density is further increased the twist angle smoothly decreases and 
asymptotes to 90° reflecting the fact that each of the two zigzag layers separately favors an anti-ferromagnetic 
orientation of neighboring instantons. 

This is the picture of the phase diagram at the level when only the straight and zigzag configurations 
are taken into account. At the next step three and more layer configurations must be considered similarly 
to the analysis performed in section [5] It may happen that the transition to three layers occurs prior to 
the change-of-twist transition discovered above. It is also natural to expect that multidimensional (D > 1) 
instanton arrays carry a non-abelian configuration of twists. Abelian configurations are always preferred 
by the straight and zigzag configurations, so it would be interesting to see when exactly the twists start 
respecting the increased dimensionality of the instanton chain. 

7 Summary 

In this work we have modeled baryonic crystals of large N c nuclear matter phase via instantons of the 
effective flavor gauge field theory on probe D8 anti-D8 branes in the geometry of near extremal D4-branes. 
The effective theory was derived as a natural generalization of the Sakai-Sugimoto [7] setup and must be 
applicable to any holographic model of quenched flavor, which enjoys the limit of large 't Hooft coupling 
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A and admits instanton description of baryons. The holographic nature of the model immediately implies 
that the most general baryon configurations are not the usual 3D crystals, but rather effectively 3+1D 
structures (including the holographic dimension). We have qualitatively analyzed the phase space of such 
structures. Since working with full-fledged 3D arrays of instantons is generally a challenging task, where 
even the pertinent flat space solutions are not known, we have resorted to several toy-models: point-like 
instantons and ID lattices. Despite the simplification we believe that these models grasp the main features 
of the 3D instanton lattices. 

We first applied the point-instanton approximation to a fD chain. To stabilize the chain in ID we 
introduced an ad hoc curvature M' in transverse spatial dimensions f|45|) . Assuming M' 3> M (curvature 
in the holographic direction) we analyzed the free energy of the zigzag configuration of point charges with 
amplitude e (figure mi)) , which is the most apparent instability mode of the straight chain (e = 0). We found 
that below critical density p c = 2 ■ 3 1 1 A M yf\ / it the straight chain is the preferable configuration, while for 
p > p c a zigzag with a density dependent amplitude (figure If 1[) has lower energy. The transition between 
the straight chain and the zigzag is second order. Considering other configurations of point charges with 
up to four layers in the holographic dimension lead to the phase diagram on figure [TBI The phase diagram 
demonstrates that in general larger densities favor configurations with more layers, although precise sequence 
of transitions, e.g. (|58|) . is not necessarily a consecutive increase of number of layers. 

After the f D example we considered a simple cubic (sc) lattice of point-like instantons. For the sc lattice 
the leading instability mode must correspond to a splitting into two fee layers. This is the 3D analog of 
the zigzag configuration. We found that for the densities p > 3.04(My / A) 3 2-layer configuration is favored 
over the original 1-layer sc lattice. Similarly bec lattice must first split into 2 sc layers, while for the fee 
lattice the double-layer splitting is probably a second order transition, which breaks cubic symmetries. The 
quantitative analysis of the 3D lattices of point instantons beyond two layers is similar to that of the ID 
chains, albeit technically more challenging. 

Next we analyzed the ID chain of full-fledged instantons. Flat space instanton solution for a periodic 
fD array of instantons was known before from the work of Kraan and van Baal 60 . We found that in the 
holographic baryon model the minimum energy of the straight chain is attained for the anti-ferromagnetic 
(twist angle <f> = n) orientation of instantons. To investigate the stability of the straight chain we constructed 
a new flat space instanton solution corresponding to the zigzag of finite size instantons. The derivation was 
based on the ADHM construction and an analog of Nahm transform similarly to the derivation of Kraan and 
van Baal. Note that to have a clear interpretation of the ADHM data as instanton positions the instanton 
separation must be larger than their size. This was achieved by taking M' ^> M - the same condition 
as for keeping the fD chain stable in spatial dimension. The analysis of the free energy of the zigzag 
instanton solution lead to the phase diagram in figure [IBJ At low densities the instanton chain is straight. 
For p > f .25M\/A (5 1 / 4 times larger than the critical density p c of point instanton lattice) the zigzag with 
e = e(p) is preferable. The transition to the zigzag is second order as in the case of point charges. As for the 
instanton orientation it is anti- ferromagnetic immediately after the phase transition, but for p ~ lASM^/X 
the twist angle changes from <^ = 7rto^><7rina first order phase transition. For larger densities </> is a 
decreasing function of the density, which asymptotes to <ft — 7r /2- To see the orientation flip we constructed a 
zigzag instanton solution for general (non-abelian) twist in the first non-trivial order in a perturbative small 
instanton expansion a <C d. Remarkably in the case of the zigzag the abelian twist are always favored over 
the non-abelian ones. 
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Analysis of the 3D instanton lattices is more challenging and was not attempted in this work. However 
the study of the above toy-models provides us with a qualitative picture of the phase space of the 3D lattices 
of holographic baryons. Namely, when the lattice is squeezed it gets stratified into 3D layers stacked one 
upon each other in the holographic dimension. Once the density is varied gradually from zero until some 
large values at which it still makes sense to talk about individual baryons the lattice undergoes a sequence 
of phase transitions with increasing number of layers. Such a behavior can already be seen at the level of 
point-like instantons. The phase space of the finite size instantons is reacher. In particular the energy of 
the instanton configuration depends on the mutual orientation of instantons. Apart from the layer splitting 
transitions one will have transitions with a change of the orientation. Another consequence of the finite 
instanton size is an effectively weaker interaction as compared to the point charges due to the screening in 
the instanton cores. 

The interpretation of the above results may be an interesting physical problem. The main idea have 
been already suggested by Rozali et al. in the case of uniform instanton distribution [18} . The stratification 
of the lattice in the fourth dimension must be a holographic analog of the emerging Fermi sea. Indeed the 
holographic coordinate plays the role of an energy scale. The energy scale associated with finite density of 
baryons must be the chemical potential. So baryons piling-up on top of each other and acquiring higher 
energy scale are filling the Fermi sea of available quantum states. In this case the Fermi sea cannot be 
baryonic: the latter are classical in this picture and carry no spin, but rather a Fermi sea of quarks. In 
a realistic setup (with no additional scale M') the critical density of the stratifying phase transitions and 
the instanton size will be defined by the same scale MvA. Therefore the transitions will happen when 
the baryons start overlapping. Overlapping baryons will loose their identities: the constituent quarks will 
no longer know which baryons they belong to. In the large N c QCD this kind of matter is believed to be 
quarkyonic, that is to have a quark Fermi sea but baryonic Fermi surface. The above phase transitions must 
then signal the onset of the quarkyonic phase 

For large enough densities the lattice will have many layers in the holographic dimension, i.e. it will 
become mesoscopic. In this regime baryons will be well overlapping, quarks will form a kind of a liquid and 
the uniform baryon density will be a good description. In such a limit the baryon popcorn picture should 
match the story of the expanding Fermi sea of Rozali et al. 

Let us finish the discussion by summarizing interesting open questions. 

• The analysis of this paper has not been generalized to the case of 3D baryon lattices. It will be 
interesting to construct new 3D instanton solutions and study the corresponding phase space. 

• Finite density baryon lattices were previously extensively studied within the framework of Skyrme 
model. Since skyrmions have a natural connection to holographic baryons there should be a way to 
reproduce the results from holographic baryon point of view. In particular one can try to understand 
the realization of the skyrmion-half-skyrmion phase transition along the lines suggested by Rho et 
al. [20]. 

• Finite temperature nuclear matter transition was studied in |17j and |18j in the approximation of 
a uniform instanton density. It would be interesting to incorporate finite temperature in the above 
analysis. 

11 Similar conclusion has been made in the case of skyrmions in [651 - The transition from the skyrmion to half-skyrmion matter 
is believed to signify the onset of the quarkyonic phase as well. Following the ideas of |20| we expect that the holographic popcorn 
transition is related to the half-skyrmion transition in Skyrme lattices. 
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• The generalization of the Sakai-Sugimoto model used in this paper is applicable for a wide class of 
holographic models of baryons. Since baryon interactions are repulsive in the Sakai-Sugimoto model 
it is desirable to make a similar analysis with attractive baryons, e.g. in the Klebanov-Strassler 
geometry [51]. It is also interesting to discuss baryons in the setup of non-critical models with A ~ 1. 
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A Straight chain. Analytical analysis at small and large densities 



In section 16.1.11 we described the flat space instanton solution that corresponds to a periodic straight chain 
of instantons. The flat space solution is the — order approximation of the holographic baryon lattice in the 
large A limit. In section [5.1.21 we used it to compute the leading 1/A corrections to the baryon free energy. 
Provided that we have a solution of the determinant of matrix L ()81|) it is straightforward to compute the 
instanton density using (|68]l and the solution of the abelian field Aq (|87|) . Formulae obtained this way are 
too bulky to present in the paper. In fact, the explicit expressions are not required for the evaluation of the 
density moments (|53")) in the non-abelian part of the energy. In the case of the abelian Coulomb energy (155)) 
however, one still needs to take an integral of a bulky expression. In this work we evaluate the integral 
analytically in the limits of small a«d and large d instantons. 

For large lattice spacing d 3> a we can expand the expression for det(L) in powers of a/d. In this case 
the instantons in the chain are approximately spherically symmetric in 4D, so at distances R ~ a <C d from 
the center we have 



logdet(L) « const + log(a 2 +i? 2 ) 



a 2 - R 2 cos(2t/>) 
~M 2 
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+ 0[a 4 /d 4 ) , (165) 



where R is the 4D radius and ip is the 4D 'latitude' angle (i. e. 
instanton density becomes 
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and consequently the Coulomb energy per instanton 

Eq 



evaluates to 
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The first term in this expression is the 4D Coulomb self-energy of the instanton, which diverges in the limit 
of zero size. The second, a-independent term is a modification of point charge energy (|51|) to include the 
dependence on the twist angle. In the d 3> a regime the Coulomb energy apparently prefers the instantons 
to be oriented in an anti-parallel way <fi = n, while the non-abelian energy (|86[) favors <f> = 0. Combining 
the two energies together in section 16.1.21 we found that the Coulomb energy wins and the equilibrium twist 
angle is (f> = it. 

In the opposite limit d < a of densely packed overlapping instantons, the instanton density becomes 
independent of the £4 coordinate, while in the 3 transverse dimensions it becomes concentrated at two 
widely separated points 



,2 



= ( O ' o 'ib (2 ^ 0) ) and * (2) = (168) 

Indeed expanding logdet(L) in powers of d/a gives 

, , 6 2-7T — 6 , (n + r 2 + (ira 2 /d)) 2 

log det(L) « 5 x ri + — ^ x r 2 + log ^ '-U- , 169 

d d or\r2 

where r x — \x — X^'\ and r 2 = \x — cf. equations (|52")l . Taking the double Laplacian one can obtain 

the corresponding 3D instanton density as 

I 3D (x) = dxI 4D (x) a ±xS^(x-X^) + ^-±x8^{x-X^) + 0{e-^ r ' d ). (170) 

Z7T Z7T 

Consequently, the abelian gauge field has 3D Coulomb form 

io = ^ + • (171) 
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Evaluating again the Coulomb energy using (|88|) gives 
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where pi and p 2 are effective radii of the instanton density concentrations that appear 5-like in equation (|170p . 
Those 5-functions are artefacts of the approximation (|169j) that becomes inaccurate near X^> and X^ 2 \ A 
better approximation for r\ ~ d <C b is given by 
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where a; = </>ri j d. The Coulomb energy of this charge density interacting with itself (but not with the other 
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charge at the X^) corresponds to 
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where Ci and C 2 are numerical values of the dimensionless integrals: 
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Likewise near X^ 2 \ 
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Plugging these charge radii into equation (jl 72[) , we write the Coulomb energy of large instantons as 
N c (C^ 2 +3(^-0) 2 ) , (C 2 + 2)^ 2 + (C 2 -2)(7r-0) 2 
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Two terms in the above expression favor different twist angles, but the first term is apparently stronger in 
the d <^ a regime. Combining this Coulomb energy with the non-abelian energy (|86j) one can find that 
(f) = 7r is a preferred angle for the overlapping instantons as well. Below in the appendix [B] we argue that 
anti-parallel orientation is favored for any values of a and d. 



B Straight chain. Equilibrium configuration 

In section 16.1.21 and appendix [X] we computed the free energy of the straight twisted instanton chain in the 
limits of small and large instantons. Here we are going to argue that the equilibrium value of the twist angle 
is always 4> = n, no matter how is the instanton size compared to the lattice spacing. First, thanks to the 
invariance of the gauge fields under Z 2 C SU(2), the instanton chain has an exact <j> — > 2tt — 4>, x% —> —X3 
symmetry. Indeed the combination of the latter transformations leave the ADHM data Y invariant up to 
an overall sign, which does not affect the gauge field. In principle, this symmetry could be spontaneously 
broken, but it does not happen for either large or small d, and it is extremely unlikely that such spontaneous 
breakdown would happen only at the intermediate d ~ ciq. Consequently, for any given d and a, the minimum 
of the net energy with respect to the twist angle <f> should lie at one of the fixed points of the transformation 
4> — > 2tt — 4>, i.e. at (f> — tt or (f> = 0. 

We are going to show that (j> = is a local maximum of the net energy Eq + £^na for the straight chain 
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instanton solution derived from the following expression for determinant of the matrix L: 
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The non-abelian energy -En A is given by equation (186[) . while the Coulomb energy must be computed from 
the integral ([55]). 

Suppose the radius a lies at the minimum of the total energy (for fixed d and <p = 0). In this case 
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From this relation at </> = and equation (|182[) one derives that 
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where the second equality follows from equation (|86j) . At the same time, 
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The derivative of the net energy with respect to <p is always negative for the equilibrium value of the radius a. 
Thus, <fi = is always a local maximum of the net energy. Since we do not expect the symmetry (f> — > 2ir — <fi 
to be spontaneously broken, the only remaining choice for the minimum is to lie at <f> = tt. 

To complete this argument it remains to derive the relation (|183|) . For this let us calculate the derivative 
dEc/dcj) for = 0. The Coulomb energy (|88|) is quadratic in logdet(L), so we need both the derivative 
ddet(L)/d(j) and the determinant det(L) itself in the 4> — > limit. Note that in this limit, the determi- 
nant (|180p becomes spherically symmetric in 3D, 
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where r 



x\. Also d'Alembertians of this expression are spherically symmetric. In such a case 
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only the spherically symmetric part of the derivative d det(L) / dcf> contributes to the integral 
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Extracting the spherically symmetric part from the expression for the derivative one gets 
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where O is the azimuthal angle in 3D. The term proportional to cosG will vanish after integration over the 
spherical coordinates. Thus, we obtain the following relation between the <j> derivative of det(L) and its 
derivative with respect to instanton radius a, 



d 2 n ddet(L) 
Air deb 



0=0 



7T 2 a 4 271T 



27ra 3 <9det(L) 



(190) 



Equation (|183|) immediately follows from the above relation between the derivatives. This completes the 
argument. With the only assumption that the Zi symmetry related to the exchange of 4> ~^ 27r — cf> is not 
spontaneously broken (which is highly unlikely since it is not broken in both a <C d and a d limits) we 
conclude that at equilibrium the straight chain always favors the anti-ferromagnetic (4> — 7r) orientation of 
instantons. 



C Zigzag - Stability at high densities 

In section 16.21 we found that for large enough densities the straight chain of instantons becomes unstable 
against the zigzag deformation (figure I10|) . This result was derived in the approximation of small instanton 
size a <§; d. The latter condition can be ensured by the high anisotropy limit M' ^> M. For M' —> M the 
critical spacing of the zigzag transition d c — > a, which is outside of the small instanton approximation. The 
analysis of the stability needs to be done numerically. It is possible that the phase diagram is qualitatively 
different for M 1 ~ M. By analyzing the opposite limit of the instanton radius «>(iwe are going to show 
that as far as only the zigzag and the straight chain are compared the former has lower free energy only in 
some intermediate density regime, while the straight chain is a preferable configuration also at asymptotically 
large densities. 

In the case of the zigzag with the antiparallel orientation of instantons we derived the following solution 
for the determinant of the matrix L in section 16.21 



det(L) / nri ira 2 nri \ ( irr2 ira 2 ttt2 
cosh — — + — — smh — — x cosh — — + — — smh 



const \ d dr\ d J \ d dt2 d 

r? + r\ — (ira 2 Id) 2 , irri . irr2 2ttx4 

+ smh — — smh — — — cos — — 

2r\T2 d d d 

7ra; 4 ( ( n , 7TT\ ira 2 irri \ ( irr 2 to 2 7rr 2 

sm v x cos — — x 2 cosh — — + — — smh — — — 2 cosh — — + — — smh 



d \\ d dr\ d J \ d dr2 d 

+ sin v X — 1 + cosh — — cosh — smh — — smh — — , (191) 

d d 2r\T2 d d J 
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where 



1,2 



(x 3 T 5& e ) 2 , 



V4e 2 + (7ro 2 /d) 2 , 
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ira 2 /d b e 

In the regime of densely overlapping instantons (d <C a) equation (|191[) yields 



(192) 
(193) 

(194) 



logdet(i) 
Hence the instanton density 
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Similar to the straight chain with c6 = it, half of the zigzag's net instanton number is concentrated at 
Xi 1 ) = (0,0,+6 e /2), and the other half at = (0,0, -b € /2), where the separation b e is given by ([153]) , 
Consequently, the abelian Coulomb energy of the zigzag has a 3D form 



Eq 



ttN c 
8XMd 



2 x 



2p 



0(djb\ 



(197) 



where p ~ d is the effective charge radius of the instanton concentration at X^ 1 ) or X^ 2 \ which appears 
8- like in equation (|196p . Note that both p and b e depend on the zigzag deformation e. To find p, we need to 
replace (|195j) with a better approximation for ri ~ d & e , namely 



logdet(L) = + log 
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7r 2 a 2 sinh(7rri/<i) 
Consequently, the instanton density near X^ comes out to be 
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sin 1/(2 + cos i/) cos(7TX4/d) 



16a 2 c? 2 (l + cosj/) sinh 5 (7rri/d) 



2 + cosh(2Wrf) - 5 >< tanll(7rri/d) 

3 + (7rri/cf) 2 sinh(27rri/<i) 
+ 2 X (27rn/d) 

-9 + 16(7rn/d) 2 + cosh(47rn/rf) 

+ 8[1 + (vrri/c?) 2 ] x cosh(27rri /<2) 

- 24(7rri/d) x sinh(27rri/d) 



Note the ^-dependent subleading term on the last line here. Its periodicity in X4 is 2c? rather than d, 
which is a characteristic of a zigzag rather than a straight-chain arrangement of the instantons. Indeed, the 
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instanton density near has a similar term but with the opposite phase. 

The Coulomb self-energy of the instanton density (|199|) near X^ 1 ' corresponds in 3D terms to a charge 
radius 

1 ird C 2 d 
2p 8d Aita' 



x (2 — cos v) . 



(200) 



where C\ and C 2 were introduced in equations (|1 T6[) and (| 1 7T[) . Overall the instanton density is concentrated 
in two tubes parallel to the X4 direction. The tubes are separated by a distance b e along the X3 direction. The 
3D radius of the tubes scales as p ~ d -C 6 e . The density is independent of £4 in the leading approximation, 
but it has a subleading modulation with the period 2d which goes in opposite phase for the two tubes. 
The net Coulomb energy (|197p of the whole zigzag is 
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Combining it with the non-abelian energy (|113p . we find the net energy cost of the zigzag deformation 
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The coefficient of e 4 is always positive, but the coefficient of e 2 can be either positive or negative, depending 
on the lattice spacing d and the instanton's formal radius a, which suggests another second-order phase 
transition between a straight chain and a zigzag, this time for d <C a. Note that the coefficient of e 2 is 
positive in the limit d — > 0, that is the straight chain becomes stable again at large densities. For the 
moment let us take this transition seriously and analyze it in more detail. 

To study the zigzag phase, it is convenient to express the energy in terms of a and 



= + (2dA) 2 x e 2 
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In the zigzag phase e ^ 0, the variables p and a are independent, albeit p > a 2 , while for the straight chain 
there is a constraint p — a 2 . In terms of these variables the non-abelian energy reads 
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and hence the net energy of the zigzag is 
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Minimizing this energy with respect to the independent variables p and a, we obtain 



^/8 



[A 2 M 2 (M /2 + |M 2 )]V4 



= >/5C 2 /8 x « 1.024 a^. (207) 



Thus in the zigzag phase the formal radius a = a z of an instanton does not depend on the lattice spacing 
d. Also, despite the d <C a regime in which individual instantons merge into two continuous lines in the £4 
direction, a z is remarkably close to the radius a' (|99p of a standalone instanton! 
The zigzag amplitude e is obtained as follows 




where a z is given in equation (|207|) and 



n = V^r _ / ^_ Xfl K 8.5-xa z (for M' > M). (209) 

V2-C 2 V M' 2 + ±M 2 Af' ; y ' 

For d < D formula (I208[) breaks down because the solutions (|206p and (|207[) fail to satisfy the constraint 
p > a 2 . For smaller lattice spacings d < D, the instantons form a straight chain, e = 0, p = a 2 , and the net 
energy should be minimized subjected to this constraint. This leads to a cubic equation for a 2 , namely 

D 2 

(2-C 2 Wxa 6 + 2C 2 a 2 z xa 4 - (2 + C 2 )a 6 z = 0. (210) 
cr 

Note that the solution is smaller than a 2 for any d < D — in the high-density straight-chain phase, the 
formal instanton radius a shrinks with the lattice spacing d; asymptotically, for d -C D, 



a 6 /2 + C 2 (d 



1/3 

X ( £ ) • (211) 



a z V2-C 2 \D 

The situation is summarized on figure I17[ where we plot the dependence of the a parameter and the 
zigzag amplitude e on the lattice spacing in the d -C a regime. Note that this whole picture — the curves 
a(d) and e(d), the transition point D, and even the existence of the transition from zigzag for d > D to a 
stable straight chain for d < D — is based on the Coulomb energy calculated in the d <C a approximation. 
Consistency of this analysis requires D <C a z , which in light of equation (|209[) means M' » 8.5M. Again, 
as in the case of small instantons, consistency requires M' M. 

Thus as far as just the straight chain and the zigzag with (f> = tt are concerned one can have the following 
possibilities. For the highly anisotropic setups with M' 3> M, the instantons form a stable straight chain 
when the lattice spacing is either very large or very small. At the intermediate lattice spacings, the straight 
ID lattice is unstable and the instantons form a zigzag. For less anisotropic setups, the zigzag region 
shrinks. For M' —> M, both transition points between straight-chain and zigzag phases move into the d ~ a 
region where our analytical calculations become unreliable. These observations are summarized in figure IT8l 
The big gray blob in the center indicates the region, where we do not have analytical control. Logically 
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Figure 17: Equilibrium instanton radius a and the zigzag amplitude e in the regime of large density a 3> d 
in the case we allow only the straight chain and zigzag configurations with 4> = n in the phase diagram. 

there are two possibilities for what is happening in this region. Perhaps there is similar phase structure 
for all M' > M: stable straight chain for either large or small lattice spacing, but a zigzag instability for 
intermediate d. Alternatively, for small enough anisotropy M'/M, the zigzag phase may disappear altogether, 
and the straight chain is stable for all lattice spacings. 

In principle one can be concerned about the meaning of the straight chain configuration at asymptotically 
large densities. In reality however we cannot just restrict to the zigzag and straight chain and must consider 
other relevant instanton configurations and compare their free energy. In particular favored configurations 
at large densities must be multi-layer 2D instanton lattices. In fact as we show in section 16.31 assuming 
M' ^> M, at high density the more stable configuration is the zigzag with tj> < tt rather than the straight 
chain, see also figure [T5] Notice that the transition to the phase with the twist angle 4> < tt occurs at critical 
spacing d ~ l/(y/XM), while the scale for the straight chain transition considered in this section is 



One can thus be worried that for M 1 ~ M, when both critical densities are of the same order, there is a 
transition from the zigzag back to the straight chain. Such a transition seems to be counterintuitive and 
would contradict the physical interpretation of the popcorn phase transitions given for example in section [7] 
Therefore it is still interesting to make a numerical analysis of the M' ~ M case and show that at high 
densities the straight chain is always dominated by multi-layer configurations. 

D Zigzag. Details of calculations for general twist 

In section HTSI we discussed the zigzag configuration with a general (non-abelian) orientation twist (|135l) , (|136|) 
between adjacent instantons. For the general twist it is also a straightforward procedure to derive the full 
instanton solution, or at least the matrix L, which is related to the topological density ([55]). However, unlike 
for the abelian limit, there is normally no way to present the results of such a calculation in a sensible 
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Figure 18: Phase diagram at different anisotropy ratios M'/M in the case we ignore multi-layer instanton 
configurations other than the zigzag with <f> = it. 

compact manner. On the other hand we have to resort to approximations anyway when evaluating the 
Coulomb energy integral (f88|) . Therefore for our purposes it is enough to calculate the instanton solutions 
up to few non-trivial orders in a/d, or its inverse. Here we will explain the calculation of the solution as an 
expansion in the case of a <C d. We will find the expression for the determinant of matrix L and compute 
the Coulomb potential Ao, instanton density / and finally the net energy of the zigzag. 

The preliminaries of the calculation were explained in the beginning of section 16.3.21 In particular we 
start from writing the matrix 

L = [x 11 - Y ti {9)) 2 + T(0) (213) 



in terms of the ADHM data (|134[) . (| 138[) defined as differential operator acting on functions ip(6) defined on 
a circle. The expansion of the ADHM data in powers of a/d in section IB. 3. 21 took the following form. The T 4 
component is exact: 
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The matrix T becomes 



(214) 



T(6) 



— (l + Eicos/3) x U(p+Z) + S(0. 



(215) 



while the remaining 3 components of r M take the form 



r(0) = en 3 E3 + r x (0) + 0(a 4 /d 3 ). 



(216) 
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Here we wrote the operators in a basis of two-component wavefunctions 



_is ' si ' * W2 ™-* /2) - <2i7) 

and Si, 2, 3 are Pauli matrices acting on the two components. The first term in (|2 16[) encodes the zeroth 
order abelian zigzag solution from section 15721 The second term Ti(9) is an 0(a 2 /d) leading correction, 
which one obtains from solving equation (|146[) . The diagonal 11 and 22 matrix elements of the solution are 
piecewise-constant functions of 9, from (|147|) 

( f i)dia g W = ^ (P + (P x q)^3 sin/3) x - 0(-| < * < f )) , 

where O is the step-function, while the off-diagonal 12 and 21 matrix elements are more complicated: 
from ([Ugl) for -0/2 < (9 < 0/2 

(f x ) 12 (0) = (f,);^) 

2 

= ^-("P3n 3 + v cosh(2(e/d)0) - i(n 3 x v) sinh(2(e/d)0)) , 

where v = — — - — j-p- (— px cosh((7r — 0)e/d) + (n3 x qj_) sin/3sinh((7r — 0)e/d)V 
cosh(7re/a) V / 

while for 0/2 < 9 < n/2, or -w/2 < 9 < -0/2, equation (|149p tells that 



± 2d 



-tg 3 sin/3n 3 + v' cosh((e/d)(7r - 2|0Q) N 
± z(n 3 x v') Binh((e/d)(7r-2|fl|)). 



where v' = — ((ri 3 x p)sinh(0e/d) — sin/3cosh(0e/d) 

cosh(7re/a) V 

In our notations, p 3 and g 3 are the x 3 components of the p and q unit vectors while p^ and q^ are their 
components perpendicular to the a; 3 direction; the ± sign in equation (1140)) is the sign of 9. 

In summary we need to compute determinant of the differential operator L{9) acting on the two com- 
ponent functions * (|2T7| .The operator L is defined through (f2T3]) with T and T 4 given by (|2lg) and ([214]) 
respectively, while the T 1,2,3 components are given by the expansion (12 1 6|) with the leading order correction 
T\ specified by equations (|147l) . (|148l) and (| 149[) . Consequently we will compute L as an expansion up to 
first few non-trivial corrections: 

L{x) = L (x) + L x {x) + L 2 {x) + ••• , (218) 
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g,2 



where 

L = (x 4 -idjL) + (x-e£ 3 n 3 ) 2 + ^- x [8(6 + |) + 8(9 - £ 

Li = -2x • fi + e{£ 3 , n 3 • f 1} + yxS[ cos/3 x ^J(0 + |) + <y(fl - |) 

£2 = (fi)' - 2x-f 2 + e{£ 3 ,n 3 -f 2 }, ••• (219) 

Notice that in the leading order we have removed a/d terms from T 3 , cf. (|105[) . Also the last term in (|218p 
is only 0{a 2 ) and naively should be dropped. Careful analysis shows however that this term is required for 
the consistency of our expansion in the regions close to instanton centers. Although the expression for L 2 
depends on the subleading corrections r 2 to the ADHM matrices, only the first term, quadratic in T\, will 
be important. The leading L$ operator has determinant!"^! 

Dct(i ) = const x ( cosh(7ri?i /d) — cos(-KXi/d) + smh^Ri/d) 

sinh( (? !)i?i/d)sinh((7r - (fyRx/d)) x 
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smh(<t>R 2 /d) sinh((7r - (p)R 2 /d) , (220) 



where 



R\ = x\ + x\ + (x 3 -e) 2 and R\ = x\ + x\ + (x 3 +e) 2 , (221) 
while the effect of the subleading L\ and L 2 operators on logdet(L) may be calculated perturbatively as 

logDct(L) = logDet(L ) + Tr (Lq x L x ) + Tr (Lq x L 2 ) - \ Tr (L„ X L X L Q ~ X L 1 ) + ••• 

+ 7T/2 +7T/2 

= logDet(Lo) + J — tr(c (0,0)^1 (#)) + I — tx(G (6,6)L 2 (6)) 

-7T/2 -tt/2 

tr{G (6 1 ,62)L 1 (6 2 )G (6 2 ,6 1 )L 1 (6 1 )) + ••• , (222) 

where Go(9i,9 2 ) is the Green's function of the operator L . At generic points in space this Green's function is 
rather messy, but for small instantons (of radius much smaller than the lattice spacing, a <^ d) the instanton 
density is concentrated in small volumes around the instanton centers where Go has a much simpler form. 



12 In the notations of this section, 'Tr' and 'Det' are the trace and the determinant in the Hilbert space of ^(9), while 'tr' 
and 'det' are the trace and the determinant of the 2x2 matrices only. 
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Specifically, at 0(a) 4D distance p from the center X 1 - 1 = (0, 0, (— l) n e, nd) of the instanton at position n, 

Go 02) = e - m < 01 - 02 U 91 °) for even n, 



, 5 2 / 

G o (0i,0 2 ) = e^Pi-to) f 52 ° 



for odd 



2 

where g 1 (8 1 ,d 2 ) = t (( independent of 6»i , 6> 2 )) + 0(l/d 2 ) 

p l + a 2 

//> « % t sinh((7r — 2I6*! — 6 2 \)e/d) ^, o , «k 

and 5 2 0i,0 2 = — x ^ ' 2V ' 1 + 0(a 2 d A ). 223 

2ed cosh(7re) 

Note that the gi is 0(l/a 2 ) while the g 2 is 0(l/d 2 ). Consequently, near instanton centers, both the first 
and the second orders of the formal perturbative expansion (|222[) may contribute 0(a 2 /d 2 ) terms to the 
determinant (|220[) . In terms of matrix elements of L\ and L 2l at distance p ~ a from the center of an 
even-numbered instanton we have 

i Dct ( L ) 2 r«w 2 r / 



2ed ./ 2tt 
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p 2 + a 2 



*/ ^(i 2 )ii(0) 



2 7T 
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+ 0{a 4 /d 4 ); (224) 



x / / -7^T sinh ((^ " 2|0i - e 2 \)e/d) x (iOu^i) x (LOaiCfe) 



near the center of an odd- numbered instanton we have a similar expression with the 22 and 11 matrix 
elements exchanging their roles. 

The diagonal matrix elements of the L\ follow from those of the V\\ in light of equation (|147|) . 



(£1)11,22 = -2(xTen 3 ) • (p cos /3 T (p x q) sin/3) x f r _ e(_| <#< Z)\ , (225) 

hence 

+ 7T/2 +7T/2 

/J/0 f c\0 

^(ii)n(0) = y ^(ii) 22 (0) = 0, (226) 

-tt/2 -tt/2 

which eliminates the first three terms on the right hand side of the expansion (|224|) . Moreover, in the 
remaining two terms we need only the leading 0(a i /d 2 ) terms in (£2)11 or 22 and the leading 0{a 2 ) terms 
in (£i)i2or2i- Since we are working at 0(a) distance from an instanton center, terms involving to x =p en3 
carry extra powers of a so that they may be neglected at this order of perturbation theory. In other words, 
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we may approximate x « ±en3, thus 



(227) 



regardless of the second-order ^(0) while 
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leads to 
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(232) 



where the last line changes sign between odd-numbered and even-numbered instantons. We are interested in 
the average energy per instanton, so in the first order of the perturbation theory that line does not contribute. 
Substituting the integrals (|230j) and (|232[) into equation ()224j) and making use of + p 2 = q^ + q 2 = 1, 
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we obtain 
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Somehow, the dependence of the intermediate expressions such as (|230[) or (|232[) on the unit vectors p and 
q cancels out from this formula, at least at the 0(a 2 /d 2 ) level. 

To obtain the instanton number density I(x) and the abelian Coulomb potential A (x) it produces, 
we need to add the perturbative correction (|233[) to logdet(Lo) and take d'Alembertians. Since we are 
working at distances p ~ a from an instanton center, we start by expanding the logarithm of det(Lo) from 
equation (|220[) in powers of a/d and p/d, 
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where Ax^ is the displacement away from the instanton center. Combining the two logarithms (|234[) 
and (|233p and taking the d'Alembertians, we find 
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and the instanton density itself 
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d 2 



(a 2 + p 2 ) 5 



d 2 (a 2 
0(l/d 4 ). 



p 2 ) cos 2 ,Ssinh((2</>- 7r)e/d) + sin 2 /?sinh(7re/d) 

X 



„2\5 



(tre/d) cosh(7re/d) 



(236) 
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Consequently, the Coulomb energy (per instanton) of the zigzag configuration is 



E c 



N c 



d 4 xA (x) x I(x) 



near center 



AM 



1 

5^ + 
4 sin 2 j3 x 



1 + 21 



2<j> 



5sinh(7re/d) - 4 sinh((2</» - w)e/d) 



80d 2 v v ^ 

sinh(7re/o?) — sinh((20 — 7r)e/d) 



(tre/d) cosh(7re/d) 
+ 0(a 2 /d 4 )] , 



(237) 



(ire/d) cosh(7re/d) 

or in terms of the generic formula (|125|) for the Coulomb energy of the instanton configuration in section [6.3.1l 

2cj>\ 2 5sinh(7re/d) - 4sinh((20- n)e/d) 



C(<f>,0,e/d) 



so 



1 + 21 



(ne/d) cosh(7re/d) 
2 sinh(7re/(i) — sinh((20 — Tr)e/d) 



4 sin f3 x 



{ne/d) cosh(7re/d) 



(238) 



Note that this Coulomb energy depends on the twist angles <f> and f3 but it does not depend on the directions 
of the twist axes p and q! 

Now consider the non-abelian energy per instanton 



E NA = XMN C { dx 4 { d 3 x (l + M 2 (x 3 ) 2 + M' 2 (( Xl ) 2 + (x 2 ) 2 )) x I[x). 



(239) 



Since the instanton density I{x) is a total derivative, we can take this integral by parts and obtain an exact 
formula for the non-abelian energy in terms of the T and T operators, namely 



E NA = N t 



\m(i + (M' 2 + ^M 2 ) x Tr(T) 



M 2 x Tr 



(n 3 -f) 2 ) 



M' 2 x Tr (f 2 - (n 3 • f)' 



(240) 



The T operator has a twist-independent trace 

-tr(T) = a 2 , (241) 
so our equation (|240[) agrees with the generic formula (|124l) for the non-abelian energy, provided we identify 

E 2 + ^xA + 0(a e /d 4 ) = Tr((n 3 -f) 2 ) ee J g tr ((n 3 • f(9)) 2 ) , 

£ x B + 0(a 6 /rf 4 ) = Tr (f 2 - (n 3 • f) 2 ) = / ^ tr (f 2 (0) - (n 3 • f(0)) 2 ) . (242) 

In perturbation theory, the bare zigzag T = eS 3 n 3 is in the a; 3 direction, hence in the transverse directions 

Tr(f 2 - (n 3 -f) 2 ) = ^(^-(na-fO 2 ) + 0(a 6 /d 4 ) (243) 
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and therefore, in light of equations (|147l) . (I148j) , and (|149l) , 

rf 2 



6 = ^xTr(f 2 -(n 3 -f^ 2 



a 1 



7T - 



4 

r 2 



x (p 2 ± cos 2 (3 + (px q)i sin 2 (3 

— x fv 2 smli^e/d) + v' 2 sinh(2(7r - 0)e/d 
8(7re/a) V 



4 



x (pi cos 2 j3 + (p x q)j_ sin 2 (3 



7r tanh(7re/d) ,_, 2 _ 2 



(ne/d) 



(pi + qisin 2 /3) 



7T 2 sinh((20-7r)e/d) 



pi - qisin 2 /?). (244) 



8 (7re/d) cosh(7re/d) 
In the £3 direction the perturbation theory is more involved: 

Tr((n 3 -f) 2 ) = e 2 Tr((S 3 ) 2 ) + 2eTr (s 3 (3 3 • f 1)) + Tr((n 3 -f^ 

+ 2eTr(l] 3 (n 3 -f 2 )) + 0(a 6 /d 4 ) , (245) 

where 

,2 



and Tr^^s-fi)) = J ^- tr (s 3 (m ■ f x ] 



= 0. (246) 

Hence 

^ x .4 = Tr ((n 3 • fi) 2 ) + 2eTr (l] 3 (n 3 ■ f 2 )) , (247) 

where the second term involves the second-order (in a 2 /d 2 ) perturbative correction to the T operator that 
we have not calculated yet. Fortunately, the A function of the moduli gives a subleading contribution to the 
net energy in the M' ^S> M limit, cf. equation (|130[) . so we do not have to work it out. 
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Altogether, in the M' 3> M limit the net energy per instanton conies out to be 



Enct = Eb + 80(XM 2 ^) X J '( moduli ) + 0(N c M 2 /M>), (248) 

where T = 3 - ^ x (2 - pi cos 2 /3 - (p x q)i sin 2 /3 

+ x ( 5 - 4sin 2 /3 + 2p 2 ± + 2qi sin 2 /?) 

(7re/a) V / 

+ f h / ( if "J )£ /,1 * (- 4 + 4sin ^ + ^ - ^) 
(7re/a) cosh(7re/a) V / 

80 



+ —r (XM 2 d 2 ) 2 x (7re/d) 2 . (249) 



In this appendix we presented the derivation of the net energy per instanton of the zigzag configuration 
with general SU(2) orientation twist in the limit of small instantons a <^ d and high anisotropy M' 3> M. 
We have demonstrated that the result can be presented in the general form of (|130p and specified the 
generic function J- in terms of zigzag amplitude, lattice spacing and twist parameters. In section \6. 3. 31 we 
minimize the net energy with respect to parameters at given spacing d to describe the phase space of zigzag 
configurations. 
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